Characterization of optically-selected star forming knots in (U)LIRGs 
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^. ABSTRACT 

Q , We present a comprehensive characterization of the general properties (luminosity functions, mass, 

H ' size, ages, etc) of optically selected compact stellar objects (knots) in a representative sample of 32 low-z 

C/3 , Luminous and Ultraluminouos Infrared Galaxies, (U)LIRGs. It is important to understand the formation 

^ ' and evolution of these properties in such systems, which represent the most extreme cases of starbursts 

in the low-z Universe. We have made use of high angular resolution ACS images from the Hubble 
Cp ' Space Telescope in F435W {~ B) and F814W {~ I) bands. The galaxies in the sample represent different 

interaction phases (first contact, pre-merger, merger and post-merger) and cover a wide luminosity range 
(11.46 < log (LiR /Lq) < 12.54). With a median size of 32 pc, most of the nearly 3000 knots detected 
?~^ , consists of complexes of star clusters. Some of the knots (~15%) are so blue that their colors indicate a 

young (i.e., < 30 Myr) and almost extinction-free population. There is a clear correlation of the mass of 
these blue knots with their radius, where M oc /ji 9i±o.i4^ similar to that found in complexes of clusters in 
f^ ' M51 and in giant molecular clouds. This suggests that the star formation within the knots is proportional 

to the gas density at any given radius. This relation does not depend significantly on either the infrared 
luminosity of the system or on the interaction phase. The star formation of all the knots is characterized 
by luminosity functions (LFs) of the knots with slopes close to 2. Nevertheless, we see a marginally 
significant indication that the LF evolves with the interaction process, becoming steeper from early to 
advanced merger phases. Due to size-of-sample effects we are probably sampling knots in ULIRGs 
C^ ' intrinsically more luminous (by a factor of about four) than in less luminous systems. They also have 

sizes and are likely to have masses characteristic of clumps in galaxies at z>l. Knots in post-mergers 
are on average larger (x 1.3-2), more luminous (2 mag) in the / band, and 0.5 mag redder than those 
in systems in earlier phases. Two scenarios are briefly discussed: (1) the likely presence of relatively 
high extinction in the most advanced mergers; (2) the dissolution of the less massive clusters and/or their 
coalescence into more massive, evolved superclusters. 

Subject headings: galaxies: evolution, galaxies: interactions, galaxies: photometry, galaxies: star clusters: general, 
galaxies: starburst 
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1. Introduction 

Star clusters tend to form where strong star forma- 
tion occurs and, specially, in st arbursts triggered by 
galax y interactions and mergers (ISchweizer & Seitzer 
19981) . During the last two decades young and mas- 
sive compact clusters have been extensively studied 
in different environme nts including star burst galax- 
ies (Meureretal.|[l995h. barred galaxies (Barthetal. 



1995 !). spiral disks (Larsen & Richder'l999'), in teract 



ing galaxies (.Whitmore et al ■1999.. ^OO?.; Bik et al. 
20031 IWeilbacher et alj l20od l2003i iKnierman et al. 



2003[ IPeterson et alJ I2OO9' Mullan et al."201l') and 
comp a ct groups of galaxies ( ^Igl e sias-Paramo & Vilche: 



2001 



Gallagher et al.l l2010t iKonstantopoulo s et al.l 



120101) . 

With masses of lO'* - 10^ M© and a median size 
of 3.5 pc (in terms of reft ), the light of these clusters 
follow the same luminosity function (LF) in the optical 
in a wide range of galactic environments, a power- 
law (^/r(L)dLocL "dL) with an index of a = 2. Yet, 
there is some controversy, since the range is rather 
large depending on the study (1.8 < a < 2.8). Recent 
studies suggest that the LF can be curved, showing a 
systematic decrease of the slope from roughly 1.8 at 
low luminosi ties to roughly 2.8 at high luminosities 
(IGieleslbOlOh . Such truncation can be reproduced by 
modeling LFs using an underlying cluster IMF with a 
Schechter function, shape of the mass function (MF) 
clai med to be observed in star-forming environm ents 



(e.g.. lHaasetZll2008llL^senll2009l:lGieIesll2009h . On 
the other hand, there are works that still favor the uni- 



versa l ity of the LF and the MF (e.g., ;W hitmor e et al. 



2007ll2010tlFall et al.ll2009HChandar et al.ll201 C 



These massive clusters do not form normally 
in isolation but tend to be clustered themselves 



(IZhang et al.ll2001t lLarsenll2004l) . Star-forming com- 
plexes represent the largest units of star fo rmation 
in a galaxy (see the review by Efremovl [l995). 
They are typically kpc-scale regions encompassing 
several clusters. Their properties have also been 
investigated in s ome nearby interac ting galaxies, 
such as in M5 1 JBastian et al.l l2005ah and in NGC 
4038/4039 JWhitmore et al.lll999b . inorder to under- 
stand the hierarchy of the star formation in embedded 
groupings. These works have confirmed that the 
properties of complexes of clusters resemble those 
of giant molecular clouds (GMCs) from which they 
are formed. In particular, the mass-radius relation, 
namely Mcmc °^ ^cmc reflects the state of virial 



equilibrium in these clouds dSolomon et al.l Il987h 
and holds down to the sca le of cloud clumps a few 
parsecs in radius (Williams et al.lll995h . Complexes 
of young clusters in M51 have a similar correla- 
tion, showing the im print from the parent GMC 
(IBastian et al.l l2005al) . Likewise, young massive 
clusters with masses above 10^ Mq follow the same 



relation (iKissler-Patig et al.l 120061) . However, young 
clusters with lower mass does not ( Larsen 2004i 
Bastian et al.1 l2005bt iKissler-PatiglTalJ I2OO6I) . 



A 

star formation efficiency that depen ds on the bind ing 
energy (predicted by Elmegreen & Efremovlll997h of 
the progenitor GMC could destroy such a relation 
during the formation of young clusters. Another 
possible explanation for this lack of a mass-radius 
relation in young clusters is that dynamical encounters 
between young clusters (and gas clouds) add energy 
into the forming clus ters, thereby increasing their radii 
dBastian et al.ll2005ah . 

Although much effort has been made to understand 
the star formation and evolutionary processes in 
interacting systems, little is known for Luminous 
(LIRGs) and Ultraluminous (ULIRGs) Infrared Galax- 
ies, which represent the most extreme cases of star- 
bursts and interactions in our nearby Universe. LIRGs 
and ULIRGs are objects with infrared luminosities of 
1O"L0 < Lboi ~ Lir[8 - lOOOfimE < IO'^Lq and 
lO'^L^ < l,r \8 - lOOOu m] < IO'^Lq, respectively 



-p> ^ 



(ISanders & Mirabell 119961) . The main source of this 
luminosity is thought to be the starburst activity, 
although an AGN may also be present, and even be 
the dominant energy source in a small percentage 
of the most luminous systems (Ge nzel et al. 1998( 
Farrah et al.l I2OO3I: lYuan et al.l I2OIOI) . They are rich 



in gas and dust, and more than 50% of the LIRGs 
and most ULIRGs show sig ns of being involved in a 
major interaction/merger (e.g Suraceetal.1 19981 2000t 



'Cui et al.' '2001'; 'Farrah et al.' '200 1': 'Bus house et all 



2002; Evans et al. 2002; Veilleux et al. 200^ 

(U)LIRGs are therefore natural laboratories for 
probing how star formation is affected by major rear- 
rangements in the structure and kinematics of galactic 
disks. Establishing the general properties (luminosity 
functions, mass, size, ages, etc) of the compact stellar 
objects in (U)LIRGs as a function of luminosity and 
interaction phase would provide relevant information 
in order to understand the mechanisms that govern 



'For simplicity, we will identify the infrared luminosity as 
LiR (Lo ) = log (L,r[8 - 1000//m]) . 



the star formation and evolution in these systems. 
It is not known either whether the processes that 
(U)LIRGs undergo produce stellar objects similar or 
very different in terms of size or luminosity than in 
less luminous (non-) interacting galaxies. 

In a cosmological context, (U)LIRGs become very 
relevant as tracers of the infrared phase (i.e., the dust- 
enshrouded star-forming phase), occurring during the 
early stages in the formation of massive galaxies. 
The so-called sub-millimeter galaxies, detected at 
z > 2, are similar to the ULIRGs observed in the 
local Universe, in the sense that they are merging 
systems with extremely high rates of star formation 



dChapman et al.1 120031: iFrayer et all 120031: Engel et al. 
12010). Besides, contrary to what it happens in the local 
Universe, (U)LIRGs are major c ontributors to the star 
formation rate density at z~l-2 (I Perez-Gonzalez et al. 
12005 ). Star-forming galaxies become increasingly 
irregular at higher redshifts, with a blue clumpy 
structure, asymmetry and lack of central concen- 
tratio n 
19961: 



(Abraham etal J Il996t Ivan den Bergh et al. 



Im et al. 1999). It has been claimed that 



clumps in clumpy galaxies represent star-forming 
complexes intrinsically more massive by one or two 
orders of magnitu des than simila r -size complexes 
in local galaxies i Efremov| 1 19951: lElmegreen et al. 



20091: [Gallagher et al. 2010). Although most of the 



clumpy galaxies, observed at z = 0.7-2, do not show 
signs of interactions, given the extreme star formation 
properties in local (U)LIRGs it is relevant to study 
the properties of their clumpy structure in order to 
establish the similarities to and differences from high- 
z clumpy galaxies. 

Previous studies of compact stellar structures in 
(U)LIRGs were focused only on small samples, mostly 
with low angular resolution ground- based imaging 
JSurace et al.lll998l:ISurace et al.ll2000l) . or on detailed 
multi-frequency studies of compact stellar obj ects in 



individual galaxies (Surace et al.-2000:Diaz-Santos et al 
20071) . Thus, no study has been done so far on a 



representative sample of luminous infrared galaxies 
covering the different phases of the interaction process 
as well as the entire LIRG and ULIRG luminosity 
range. 

This paper presents the first attempt at obtaining an 
homogeneous and statistically significant study of the 
photometric properties (magnitudes, colors and sizes) 
of optically-selected compact stellar objects found in 
these systems as a function of Lm, morphology (i.e., 
interaction phase) and radial distance to the nucleus 



of the galaxy. These stellar objects, some of them 
identified as YMCs candidates for the closest galaxies 
(i.e., at Dl <60 Mpc) and practically all of them 
identified as cluster complexes further away, will be 
referred to as "knots". The luminosity function of 
the knots is also evaluated and compared with that 
obtained in other studies. This study also represents 
the starting point of a more quantitative analysis based 
on the direct comparison of these measurements with 
state-of-the-art simulations of galaxy encounters with 
linear resolutions comparable to those presented here 
(e.g., simulations by Bournaud et al. 2008b). The 
results of this analysis, in particular the evolution of 
the luminosity functions of the knots, will be presented 
in a forthcoming paper (Miralles-Caballero et al. 20 11 , 
in preparation). 

The paper is organized as follows. We first describe 
the sample and classify the systems according to their 
morphology in section |2l The data are presented in 
section [3] We then explain the analysis methodology, 
including the source selection, photometric measure- 
ments, completeness tests, etc. In section|5]we present 
the main photometric results, including the properties 
of the objects measured as a function of Lir and the 
interaction phase of the system. Age and mass deriva- 
tion for the young objects is also discussed. Finally, 
the summary and conclusions are given in section |6] 
Throughout this paper we use Qa - 0.73, Q.m - 0.27 
and //() = 73 kms-^Mpc-\ 

2. Sample and properties 



2.1. Sample selection 

Galaxies were selected from t he flux-limited IRAS 
Rev isited Galaxy Sample (RBGS. ISanders et alj|2003h 
and ISanders et alj (Il988h . using the following addi- 
tional criteria: to sample the wide range of IR lumi- 
nosities in (U)LIRGs, to cover all types of nuclear 
activity, to span different phases of interaction and to 
optimize the linear scales by selecting low-z galaxies. 

For all the selected systems, archival high angu- 
lar resolution HST blue and red optical images are 
available. The / band, less affected by extinction, 
can roughly indicate the extent of the system and the 
tidal features emerging from it, whereas the B band 
can give us an idea of the young population. Two of 
the systems were not considered for this study, since 
they are so close (< 40 Mpc) that individual stars 
belonging to them start to contaminate, and selection 
effects can be encountered for these galaxies. For 
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Fig. 1. — : Distribution of the 32 systems of the sample in 
the luminosity-redshift plane. The sample consists of low- 
z (z<0. 1) systems and covers a factor of ~ 10 in infrared 
luminosity. 



instance, the faintest objects de tected in IRAS 10257- 
4338 by iKnierman et alj (l2003k with the WFPC2 can 
be individual stars (though they expect them to be 
few). Since the ACS images are more sensitive at 
similar exposure times, more individual stars can be 
detected. In addition, the galaxies in the sample were 
observed with IFS (VIMOS, INTEGRAL and PMAS), 
with the aim at identifying Tidal Dwarf Galaxy candi- 
datefl 

Our sample comprises 32 low-z systems, 20 LIRGs 
and 12 ULIRGs. Individual properties of the systems 
are presented in Table [T] It is not complete in either 
volume, flux or luminosity. However, this sample is es- 
sentially representative of the local (U)LIRG systems 
(see table [1] and Figure [TJ: (1) it covers the luminosity 
range 11.39<Lir < 12.54 and the redshift range 0.016 
< z < 0. 124 (from 65 to about 550 Mpc), with a median 
value of 0.037; (2) it spans all types of nuclear activity, 
with different excitation mechanisms such as LINER 
(i.e., shocks, strong winds, weak AGN), HII (star 
formation) and Seyfert-like (presence of an AGN); (3) 
and all the different morphologies usually identified in 
these systems are also sampled (see section 12.21 and 
Figure|2|. Therefore, the sample is appropriate for the 
purpose of this work. 



2.2. Morphology 

With a sample of 32 systems we have the oppor- 
tunity to investigate the properties of their compact 



stellar objects as a function of the morphology of 
the galaxies. Besides, setting a rough estimate of 
the timing of the merger according to its morpholog- 
ical class, we will compare the properties of these 
objects with those of compact stellar objects mea- 
sured in high resolution interacting galaxy simulations 
by iBournaud et al.l ( l2008bl) in a forthcoming paper 
(Miralles-Caballero et al. 2011, in preparation). The 
galaxies exhibit a variety of morphologies, from a 
wide pair with a projected nuclear distance of 54 kpc, 
to close nuclei with a projected radial distance of 
d<5 kpc, and to a single nucleus with a distorted stellar 
envelope. We derived a morphology class for each 
system of the sample by using the red (F814W filter 
in ACS and WFPC2 instruments) HST images. A 
m erging classification scheme similar to that given 
bv lVeilleux et al.l(l2002h has been considered. Though, 
due to our limited number of galaxies some grouping 
has been made, so as to have enough statistics in all 



phases. For instance, phases I and II in lVeilleux et al. 



(I2002h have been grouped for this study. The sample 
has been divided in four groups according to their pro- 
jected morphology and trying to match well identified 
temporal snapshots in the evolution of the interaction 
according to the models (see images in Figure|2]i: 

I-II. First Approach- In this early phase, the first 
passage has not been completed. It is prior to the first 
close passage or during that first passage. The galaxy 
disks remain relatively unperturbed and bars and tidal 
tails have not yet formed. In some cases the disks can 
overlap, but no projected morphological disruption is 
seen on either galaxy. Some examples of systems in 
this phase are IRAS 06259-4708 and IRAS 07027- 
6011. 



III. Pre-Merger. - Two identifiable nuclei, with 
very well identified tidal tails and bridges, charac- 
terizes this phase. Both disks are still recognizable. 
This occurs 200-400 Myr after the first approach, 
ac cording to the similar morphologies observed 



Bournaud et al.l (I2008H) models. Some examples 



These data will be presented in a second paper. 



representing this phase are IRAS 08572-1-3915 and 
IRAS 12112+0305. 



l y. Merger. - We ha ve only included phase IVa 
of lVeilleuxet al.l (120021) . since the systems in phase 
IVb have morphologies more similar to those in 
phase V than in IVa. In this phase both nuclei have 



Table 1 
Sample of (U)LIRGs 



IRAS 


Other name 


R. A. (J2000) 


DEC (J2000) 


Dl 


L,R 


Moi-phology 


Interaction 


Spechral Class 


Ref. 


name 




(hh:mm:ss) 


(deg:mm:ss) 


(Mpc) 






phase 


(Optical) 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


00506+7248 




00:54:04.91 


+73:05:05.5 


65.3 


11.46 


IP 


III 


HII 


[I] 


02512+1446 




02:54:01.99 


+ 14:58:26.2 


131.3 


11.63 


IP 


III 


HII 


[2] 


04315-0840 


NGC 1614 


04:33:59.99 


-08:34:42.8 


66.1 


11.67 


SN 


IV (2) 


HII 


[3] 


05189-2524 




05:21:01.42 


-25:21:48.3 


180.6 


12.15 


SN 


V(2) 


Sy2 


[4],[5] 


06076-2139 




06:09:45.98 


-21:40:25.6 


158.3 


11.65 


IP 


III (1) 


Sy/L-HII 


[6] 


06259-4708 




06:27:22.53 


-47:10:45.4 


164.2 


11.92 


MI 


I-II (1) 


HII-HII 


[7] 


07027-6011 




07:03:27.46 


-60:16:00.2 


131.8 


11.62 


IP 


I-II (0) 


Sy2 


[7] 


08355-4944 




08:37:01.98 


-49:54:27.5 


108.5 


11.58 


SN 


IV (2) 


HII 


[8] 


08520-6850 




08:52:28.86 


-69:02:00.3 


197.1 


11.81 


IP 


I-II (1) 




- 


08572+3915 




09:00:25.62 


+39:03:54.1 


250.6 


12.13 


IP 


III 


L/HII 


[9], [10] 


09022-3615 




09:04:12.26 


-36:26:58.1 


256.4 


12.30 


SN 


IV (2) 






F1003S-3338 




10:06:05.04 


-33:53:14.8 


143.8 


11.75 


DN 


IV (2) 




- 


10173+0828 




10:20:00.10 


+08:13:33.3 


209.4 


11.84 


SN 


V 




- 


12112+0305 




12:13:45.82 


+02:48:38.0 


318.3 


12.34 


IP 


III 


L-L 


[4], [5], [9] 


12116-5615 




12:14:21.77 


-56:32:27.6 


113.7 


11.59 


SN 


V(2) 


HII 


[7] 


12540+5708 


Mrk231 


12:56:13.73 


+56:52:29.6 


178.9 


12.54 


SN 


V(2) 


Syl 


[2], [5] 


13001-2339 




13:02:52.49 


-23:55:19.6 


90.7 


11.46 


DN 


IV (2) 


L 


[11] 


13428+5608 


Mrk 273 


13:44:41.93 


+55:53:12.3 


159.8 


12.14 


DN 


IV 


L-Sy2 


[2] 


13536+1836 


Mrk463 


13:56:02.80 


18:22:17.20 


215.0 


11.82 


DN 


m 


Syl-Sy2 


[12] 


14348-1447 




14:37:40.36 


-15:00:29.4 


361.6 


12.35 


DN 


m 


L 


[5], [9] 


15206+3342 




15:22:38.00 


+33:31:35.9 


559.6 


12.28 


SN 


V 


HII 


[5], [13] 


15250+3609 




15:26:59.05 


+35:58:38.0 


236.3 


12.05 


SN 


IV 


L 


[5] 


15327+2340 


Arp220 


15:34:54.63 


+23:29:40.5 


75.5 


12.16 


DN 


IV 


L 


[13] 


I6I04+5235 


NGC 6090 


16:11:40.10 


+52:27:21.5 


123.1 


11.50 


IP 


III 


HII 


[2] 


F17138-1017 




17:16:35.64 


-10:20:34.9 


72.2 


11.39 


SN 


V(2) 


HII 


[14] 


17208-0014 




17:23:19.14 


-00:17:22.5 


181.7 


12.41 


SN 


V 


L 


[5], [15] 


FI 8093-5744 


IC 4686 


18:13:38.66 


-57:43:53.7 


72.2 


11.55 


MI 


10(1) 


HII-HII 


[7] 


18329+5950 


NGC 6670 


18:33:37.06 


+59:53:19.3 


121.2 


11.63 


IP 


m 


- 




20550+1656 




20:57:23.51 


+ 17:07:34.6 


152.5 


11.92 


IP 


III 


HII 


[2] 


22491-1808 




22:51:45.79 


-17:52:22.5 


338.6 


12.19 


DN 


10(1) 


HII 


[5] 


23007+0836 




23:03:16.92 


+08:53:06.4 


67.9 


11.54 


IP 


I-II 


Syl-Sy2/L 


[2] 


23128-5919 




23:15:46.46 


-59:04:01.9 


189.6 


11.96 


DN 


III (I) 


Sy/L-HII 


[5], [9] 



Note. — Col (1): object designation in the IRAS Point and Faint Source catalogs. Col (2): other name. Cols (3) and (4): right ascension (hours, minutes 
and seconds) and declination (degrees, arcminutes, arcseconds) taken from the from the NASA Extragalactic Database (NED). Col (5): luminosity distances 
derived assuming a ACDM cosmology with Ho = 73 kms~^ M pc~^ , Qa=0.73 and £2^=0.27, and using the Eduard L. Wright Cosmology calculator, which is 
based in the prescription given by Wright 1 2006). Col (6): logarithm of the infrared luminosity, Lir - L(8 - 1000m), in units of solar bolometric luminosity, 
computed following Sanders & Mirabel 1996 (Lir = 4;r£)[^ Fir, where Fir = LB-vlO"''^ x (13.48/i2 + 5.16/25 + 2.58/60 + /ioo)[Wm--], being /i2j25,/60 and 
/lop the IRAS flux densities in Jy at 12, 25. 60 and 100 fiin). The IRAS flux densities were obtained from the IRAS Point Source & IRAS Faint Source catalogs 
(Moshir et al. 1993; Joint has Science 1994). Col (7): MI stands for Multiple Interacting galaxies, IP interacting pair, DN double nucleus and SN single nucleus. 
Col (8): moiphological classification used in this study, where I-II stands for first approach. III for pre-merger, IV for post-merger and V for relaxed system 
(see text). Between brackets the classification by Rodriguez-Zauri'n et al. 1 201 1} is given, when available. In their study, class corresponds to classes I-II, 
class 1 to class III and class 2 to classes IV-V in this study, respectively. Note a disagreement in IRAS 06259-4708 and IRAS 08520-6850 due to the shght 
different of the definitions of the early interaction phases in both studies (see text). Col (9): spectral class in the optical. A slash between two diff^erent classes 
indicates that no clear classification is known between both, whereas a dash separates the class of each nucleus in the system when known. L stands for LINER 
and Sy refers to a SeyfeiT activity. Col (10): references to the spectral cla ss- [1] Alonso-Heirero et al. 1 2009), [2] Wu et al. i 1998), [3] Kotilainen et al. i 2001 ), 
[41;^R isaliti e t alj t200@, [51iNar dini et_aUa008). [61 Arribas et all (200^, r7irKewlev e t al. 1 2001), [8] Cohen 1 19 92), [ 9] Evans et al. i2002), [1 0] Arribas et aU 
(2000) , [ll1 l'Zenner&Lenzer^h99jb ,^^1 lGarclVMariI^|20Q7^ .rl31 l^ [I41 ICorbett eFallfMBl) . [ISI IArribas & Colin jf2003l) . 




Fig. 2. — : Systems under study. False color composite image of the complete sample using F814W (red) and F435W (blue) //ST- 
AGS images. We have saturated the images below a given surface brightness, typically 5-10 times the global background deviation. 
The result of this, the diffuse red light, shows the lowest surface brightness features (tails, plums and shells). The white horizontal 
line indicates a scale of 5kpc. Blue knots along the tails are clearly visible in some systems. North points up and East to the left. 




Fig. 2. — : - Continued 




Fig. 2. — : - Continued 



apparently coalesced. We consider that the nuclei 
have coalesced as soon as their separation is less than 
1 .5 kpc, since numerical simulations have shown that 
by the time the two nuclei have reached a separation 
of < 1 kpc, the stellar system h as basically a chieved 



equilibrium (Mihos 1999; Naab et al.lbooel) . Long 



tails and shell structures are seen in this phase. Both 
disks are no longer recognizable, but a common 
internal structure has formed. This phase is seen 
500-700 Myr after the first approach, according 
to the models mentioned. IRAS 04315-0840 and 
IRAS F10038-3338 are representative galaxies of 
this phase. It has to be taken into account that we 
are considering projected distances. Hence, those 
systems with double nuclei allocated to this phase 
could well belong to the previous phase if the real 
distance of the nuclei is considerably greater 

V. Pos t-Merger. - We hav e also included phase 



IVb of Veilleux et al. 



(120021) . In this phase there 
are no prominent and bright tidal tails or bridges, 
since their surface brightness starts to fall below the 
detection limit, though some shell structures around 
the nucleus can still be detected. These isolated 
nuclei can have disturbed stellar envelopes, similar to 
those in the previous phase, indirectly indicating the 
past interaction. Typical ages of t > 1 Gyr after the 
first passage characterizes the timescale of this phase, 
according to the models mentioned. IRAS 05189- 
2524 and IRAS 12116-5615 are good examples of 
systems in this phase. 

Since this classification has been performed only 
by evaluating the apparent morphology of the sys- 
tems under study, it is always subject to improve- 
ments based on data with a higher angular resolution 
and sensitivity. Therefore, misclassification could oc- 
cur in a few systems. In addition, for certain stud- 
ies, slightly different and simpler classifications have 
been used previously (Rodriguez-Zaurinet al. 201 1| 
see also Table 1). Although they have the advantage of 
reducing subjectivity and therefore uncertainty, their 
classification is too coarse for detail comparison with 
model predictions. While the classification of a few 
objects can differ, as is the case for IRAS 06259-4807 
and IRAS 08520-6850 in common with Rodriguez- 
Zaurin's sample (the definitions of the early phases 
in both studies are slightly different), the agreement 
between both classifications is excellent. Using this 
classification, 4 systems are assigned to category I-II, 



13 to category III, 8 to category IV and 7 to category 
V (see table [B. 

3. The data 

The results presented in this paper are based on high 
angular resolution archival HST images. The available 
dataset for each galaxy at the moment of the analysis 
is presented in table |2l 

Advanced Camera for Surveys (ACS) broad-band 
images were taken for thirty systems from the Hubble 
Legacy Archive, with the filters F814W and F435W. 
The former is equivalent to the ground-based Johnson- 
Cousins /, whereas the latter differs from the ground- 
based Johnsoi>CousinsB between 7 and 12% in 
flux (ISiriannietalJl2005b . The field of view (FoV) 



is 210" X 210". Total integration times were taken 
from 720 to 870s with filter F814W and from 1260 
to 1500s with filter F435W. For IRAS 135364-1836 
and IRAS 15206-H3342, Wide Field Planetary Cam- 
era 2 (WFPC2) images were taken, with the filters 
F814W and F439W, being equival ent to ground-based 
Johns on / and B, respectively dOrigha & Leitherer 
2000). Here, the systems fall in the PC chip, giving 
a FoV of about 40" x 40". Whenever both ACS 
and WFPC2 images were available in both filters, 
ACS images were preferred for its superior sensitivity 
and larger FoV. The pixel size in both instruments is 
about 0.05". From now on, we will refer to the blue 
filter as F435W for ACS and WFPC2 for simplicity. 
Complementary Near Infrared Camera and Multi 
Object Spectrometer (NICMOS) broad-band images 
with filter F160W {~ H) were also available for twenty 
galaxies, which were helpful to locate the nuclei of the 
systems. However, they have not been included in 
the photometric study since they have small FoV (20" 
X 20"), not covering the outskirts of the systems or 
interacting tail-bridge structures. 

Images were reduced on the fly {HST pipeline), 
with the highest quality available reference files at the 
time of retrieval. The calibrated ACS F435W images 
had only few cosmic rays remaining, since 3 images 
were observed and then combined. However, all the 
ACS F814W images were severely contaminated with 
cosmic rays in a ~ 6 arcsec width band across the entire 
FoV. Two exposures were taken and when they did not 
overlap the HST pipeline was not able to discriminate 
the cosmic rays. In most cases that band went through 
part of the system, covering about 10-20% of its exten- 
sion. Cosmic ray rejection was then carried out using 



Table 2 
Available HST data: observing log details 



IRAS 


HST opti- 


PI and 


Obs. 


Exp. 


HST NIR 


PI and 


Obs. 


Exp. 


name 


cal data 


proposal ID 


date 


time (s) 


data 


proposal ID 


date 


time (s) 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(V) 


(8) 


(9) 


00506+7248 


ACS/WFC 


Evans, 10592 


04/09/05 


800, 1500 


NICMOS 


Alonso-Herrero, 10169 


09/10/04 


311 


02512+1446 


ACSAVFC 


Evans, 10592 


22/07/06 


720, 1260 


- 


- 


- 


- 


04315-0840 


ACSAVFC 


Evans, 10592 


14/08/06 


720, 1260 


NICMOS 


Rieke, 7218 


07/02/98 


191 


05189-2524 


ACSAVFC 


Evans, 10592 


13/08/06 


730, 1275 


NICMOS 


Veilleux, 9875 


28/08/04 


2560 


06076-2139 


ACSAVFC 


Evans, 10592 


14/11/05 


720, 1260 


- 


- 


- 


- 


06259-4708 


ACSAVFC 


Evans, 10592 


19/04/06 


780, 1350 


- 


- 


- 


- 


07027-6011 


ACSAVFC 


Evans, 10592 


19/09/05 


830, 1425 


- 


- 


- 


- 


08355-4944 


ACSAVFC 


Evans, 10592 


08/09/05 


780, 1350 


- 


- 


- 


- 


08520-6850 


ACSAVFC 


Evans, 10592 


15/04/06 


870, 1485 


- 


- 


- 


- 


08572+3915 


ACSAVFC 


Evans, 10592 


09/12/05 


750, 1305 


NICMOS 


Maiolino, 9726 


17/03/04 


160 


09022-3615 


ACSAVFC 


Evans, 10592 


07/06/06 


750, 1305 


- 


- 


- 


- 


F10038-3338 


ACSAVFC 


Evans, 10592 


06/11/05 


740, 1290 


- 


- 


- 


- 


10173+0828 


ACSAVFC 


Evans, 10592 


01/12/05 


720, 1260 


- 


- 


- 


- 


12112+0305 


ACSAVFC 


Evans, 10592 


20/02/06 


720, 1260 


NICMOS 


Scoville, 7219 


15/11/97 


192 


12116-5615 


ACSAVFC 


Evans, 10592 


07/09/05 


830, 1425 


NICMOS 


Surace, 11235 


01/07/07 


2495 


12540+5708 


ACSAVFC 


Evans, 10592 


1 1/05/06 


830, 1425 


NICMOS 


Veilleux, 9875 


09/09/04 


2560 


13001-2339 


ACSAVFC 


Evans, 10592 


26/02/06 


720, 1260 


- 


- 


- 


- 


13428+5608 


ACSAVFC 


Evans, 10592 


17/1 1/05 


820, 1425 


NICMOS 


Maiolino, 9726 


1 1/05/04 


599 


13536+1836 


WFPC2/PC 


Sanders, 5982 


27/1 1/95 


360, 1030 


NICMOS 


Low, 7213 


27/12/97 


480 


14348-1447 


ACSAVFC 


Evans, 10592 


06/04/06 


720, 1260 


NICMOS 


Scoville, 7219 


31/12/97 


480 


15206+3342 


WFPC2/PC 


Sanders, 5982 


05/09/95 


343, 750 


- 


- 


- 


- 


15250+3609 


ACSAVFC 


Evans, 10592 


25/01/06 


750, 1305 


NICMOS 


Scoville, 7219 


20/1 1/97 


223 


15327+2340 


ACSAVFC 


Evans, 10592 


06/01/06 


720, 1260 


NICMOS 


Maiolino, 9726 


10/01/04 


599 


16104+5235 


ACSAVFC 


Evans, 10592 


18/09/05 


800, 1380 


NICMOS 


Maiolino, 9726 


01/12/03 


599 


F17138-1017 


ACSAVFC 


Evans, 10592 


31/03/06 


720, 1260 


NICMOS 


Alonso-Herrero, 10169 


23/09/04 


240 


17208-0014 


ACSAVFC 


Evans, 10592 


04/04/06 


720, 1260 


NICMOS 


Scoville, 7219 


26/10/97 


223 


F 18093-5744 


ACSAVFC 


Evans, 10592 


08/04/06 


830, 1425 


NICMOS 


Alonso-Herrero, 10169 


26/09/04 


240 


18329+5950 


ACSAVFC 


Evans, 10592 


31/10/05 


830, 1425 


NICMOS 


Surace, 11235 


07/03/09 


2495 


20550+1656 


ACSAVFC 


Evans, 10592 


15/04/06 


720, 1260 


- 


- 


- 


- 


22491-1808 


ACSAVFC 


Evans, 10592 


04/05/06 


720, 1260 


NICMOS 


Scoville, 7219 


21/11/97 


480 


23007+0836 


ACSAVFC 


Evans, 10592 


12/06/06 


720, 1260 


NICMOS 


Scoville, 7219 


10/11/97 


351 


23128-5919 


ACSAVFC 


Evans, 10592 


22/03/06 


830, 1425 


NICMOS 


Maiolino, 9726 


13/10/03 


599 



Note. — Col (1): object designation in the IRAS Point and Faint Source catalogs. Col (2): optical instruments used on board HST . WFC 
refers to the Wide Field Channel of the ACS camera. Sources with the instrument WFPC2 were focused on the PC chip. Col (3): last name 
of the principal investigator, followed by the proposal ID of the optical observation. Col (4): day when the optical observation was taken. Col 
(5) total exposure times for each galaxy, with filter F814W (left) and filter F435W (right). For observations with the WFPC2 the blue filter 
corresponds to F439W. Col (6): NIR instruments used on board HST to observe the retrieved F160W images. For all these observations the 
camera used was NIC2. Col (7): same as (3) for the infrared data. Col (8): same as (4) for the infrared data. Col (9): same as (5) for the 
infrared data, with filter F160W. 
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the IRAF task credit with the aid of the blue images 
since they were barely contaminated. Most of the 
systems have observations that cover the entire FoV 
with only this instrument and the two filters, so the aid 
of multi-band images was limited to few systems. 

4. Analysis techniques 

4.1. Source detection and photometry 

Compact regions with high surface brightness 
(knots hereafter) that are usually bluer than the 
underlying galaxy have been identified for each 
galaxy. They have been detected inside a box that 
encloses each system including the diffuse emis- 
sion observed in the red band. They appear sim- 
ilar to the bright blue knot s found in interactin g 
galaxies like i n NCG 7 252 ('Whitmore et al.' '1993), 
NGC 403 8/39 JWhitmore & Schweizer. .19951 some 
ULIRGs r 



Lee et al 



Surace et al J 1 1 998h M5 1 dBik et al. 12003 . 
2005|) and Arp284 jPeterson et al.l UOO^ 



These knots were identified on the basis of having 
a flux above the 5cr detection level of the local 
background. Due to the irregular and distorted 
morphologies for most of the systems, we performed 
a two step analysis to identify the compact regions: 
we first smoothed the image, averaging in a box of 
30x30 pixels, and subtracted this smoothed image 
from the original one in order to make a first rough 
local bac kground s ubtraction; we then ran SExtrac- 
tor CBertin & Arnouts. .1996.) on the resulting image 
to detect the knots. Due to the steepness of the light 
profile of the local background in the inner regions 
and its larger deviation with respect to outer regions, 
an appropriate background subtraction was not always 
possible and some inner knots were never detected. 
Thus, visual inspection was needed afterward to 
include some inner knots and also to ehminate some 
spurious detections. 

Once the knots were detected, all photometric 
measurements were done using t he tasks PHO T and 
POLYPHOT within the IRAF ( JTodvl 1 1993b envi- 
ronment. We first identified the point-like sources 
against the extended ones, since the photometry of 
the former requires aperture corrections. We fitted 
the knots with a Moffat profile, and those with a 
FWHM < 1.5xFHWM the profile of the stai's in the 
field were considered as point-like knots. Note that 
although they actually consist of point-like and slightly 
resolved objects, they are considered as point-like 
objects here just to assign an aperture radius, since 



both require aperture corrections when the photometry 
is done. Once their size is derived via psf fitting (see 
section 14.21 ) we differentiate between resolved and 
unresolved knots. 

Aperture photometry with a 3 pixel aperture radius 
(0.15") was performed for the point-like knots. For 
the resolved knots, we used aperture radii between 
5 and 7 pixels (0.25-0.35"), depending on when the 
radial profile reached background levels. The 7-pixel 
apertures were large enough to include practically all 
the light from the knot but small enough to avoid 
confusion with other sources and to avoid adding more 
noise to the measurements. In some cases, the irregular 
shape of the knots demanded the use of polygonal 
apertures. In some cases, especially for knots be- 
longing to the most distant galaxies, the polygonal 
apertures were larger than 7 pixels. Estimates of the 
underlying background flux were made by using the 
mode of flux values measured in an 8-13 pixel annulus 
centered on the computed centroid with POLYPHOT, 
and eliminating the 10% of the extreme values at each 
side of the flux background distribution (e.g., if there is 
a nearby knot it adds spurious flux at the bright end of 
the distribution). Most of the knots are apart from each 
other by more than 7 pixels (0.35"), hence confusion 
does not seem to be a serious problem, though in 
the inner regions of few systems (IRAS 04315-0840, 
IRAS 15206-H3342 and IRAS 20550-H1656) there is 
some overcrowding. 

Photometric ca l ibratio ns were performed follow- 
ing Sirianni et alj (20051) for the ACS images and 
Holtzman et all d 1995b for the WFPC2 images (here- 
after calibration reference papers). Magnitudes were 
then computed in VEGA system. Correction for 
charge transfer efficiency (CTE) was needed in the 
WFPC2 photometry but was not applied to the ACS 
photometry, since tests in some images shown that 
the effect was less than 1% for the faintest regions 
identified. Aperture corrections were considered for 
the point-like knots. To that end, we performed the 
photometry for isolated stars in the field (from a few 
to about a hundred, depending on the image) in the 
same way as we did for the point-hke knots, then a 
second time but using a 10 pixel radius. The ratio of 
both values gives us an aperture correction up to 10 
pixels. From 10 pixels to an infinite aperture radius 
the values in the calibration reference papers were 
taken, since the PSF profile from 10 pixels on is stable 
enough regardless of the observing conditions. Typical 
aperture corrections were about 0.3 and 0.2 mag for 
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the F814W and F435W filters, respectively. Using 
the deviation of the different aperture corrections 
computed with stars in the field we estimated the 
uncertainty due to the aperture corrections to be of 
the order of 0.05-0.1 mag in both filters. Typical 
photometric uncertainties for all the knots lie between 
0.05 and 0. 15 mag, depending on the brightness of the 
knot. According to the completeness test performed 
in section |431 the systematics affect both filters nearly 
equally, thus the typical uncertainty in color is between 
0.10 and 0.20 mag. 

We corrected all the magnitudes for reddening 
due to our galaxy taking the values directly from 
the NASA/IPAC Extragalactic Database (NED), com- 
puted following Schlegel et al. (1998) . Apart from 
this, no internal reddening corrections have been 
applied to the magnitudes reported in this paper. 

Initially, we probably detected old globular clus- 
ters, super star clusters, star complexes, Hii regions, 
foreground stars, background objects, etc. Some of 
them were rejected after the photometry was made. 
Foreground stars were mostly easily identified by their 
high brightness and red colors and their status as point- 
like objects. The fields in systems at low galactic 
latitudes (IRAS 08355-4944, IRAS 09022-3615 and 
IRAS 12116-5615) are so crowded with foreground 
stars that some may still contaminate the photometric 
sample knots after rejection, especially if a faint fore- 
ground star lies within the inner regions of the galaxy. 
No cluster/knot with color Mf4-i5w-Mf^i4w ^ 2.2 fits 
in stellar population models for an extinction-free star- 
burst at redshift z~0. Since we do not expect much 
extinction in the outer parts of the galaxies (usually 
where we do not see difiiise emission from the galaxy 
in the red band), all the red knots outside the diffuse 
emission of the systems are likely to be background 
objects and were also rejected. After these rejections a 
total of 296 1 knots were considered under study. 

The nuclei of the systems were identified using 
the H band images and with the aid of other works 



that have detected them in radio (e.g., iDinshaw et al- 
ii 999i) . The identification of the true nucleus is of 
relevance in the derivation of the projected distances 
of the knots with respect to the closest nucleus, to 
properly study the spatial distribution of these knots. 



4.2. Determination of sizes 

4.2.1. Effective radius 

The half-light radius of the knots was derived by 
fitting radial Gaussians, technique that ha s been suc- 
cessfully used in other works (e.g.. lSurace et a l. 1998t 



Whitmore et al.l 119931) . The measured radial profile 



of a knot corresponds to the convolved profile of the 
PSF with the intrinsic profile of the knot. We then 
fitted 2D-Gaussians to each knot and to foreground 
stars in the image (to compute the PSF profile). We 
can estimate the cr for the intrinsic profile of the knot 



as cr- Jcrlj - o'p^f, where cr„ corresponds to the cr of 



the measured profile (average of cr in the x and y direc- 
tion) and crpsF to that of the PSF profile. An effective 
radius, r^ff, was derived as the half value of the FWHM 
derived from the computed cr (FWHM=2.354cr). 

Note that reft defined in this study does not strictly 
correspond to the half-light radius, since the light pro- 
file of the globular clusters and young massive clusters 
detected in the Milky way and in other galaxies is not 
necessarily Gaussian. However, given the fact that the 
distance of 80% of the systems is larger than 100 Mpc 
we do not expect to measure resolved sizes as small 
as the typic al size of young massiv e clusters (typically, 
r^ff < 20 pc dWhitmore et al.ll999h . but associations or 
cluster complexes. Sizes of globular clusters in the 
Milky Way are even smaller Therefore, it is no t worth 
using in this study more robust methods (e.g.,|Larsen 



ng] 



1999) that better characterize the light profile of the 
source, and where the derived r^ff is closer to the half- 
light radius . The main interest here is the relative 
change in the sizes of the knots among the diff'erent 
systems. 

The rms scatter around the 1 : 1 line in Fig|3]is about 
30%, from 3 pc to 90 pc for knots with sizes of 10 
and 300 pc, respectively. This should give a reason- 
able estimate of the uncertainty in the measurement 
of Teff. Differential extinction, which is wavelength 
dependent, is also responsible for part of the scat- 
ter. Then, under similar local background conditions, 
knots will look a bit more extended when observed 
with the F814W filter. Another source of uncertainty 
comes from the young population. In fact, faint young 
blue knots easily detectable in the F435W image but 
with low S/N ratio in the F814W image will look 
more extended in the blue ones. In any case, this 
measurement is unlikely to be affected by the nebular 
emission-lines (i.e., Hq,, Hp or oxygen lines), since 
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Fig. 3. — : Relation between effective radius for all the knots 
measured in images with filter F814W and F435W. Both 
axes are in log scale. The line indicates a ratio of unity. 



in the vast majority of the systems these lines do not 
significantly contaminate the HST filters used in this 
study. 

In order to estimate the minimum size we can 
achieve by using this technique we estimated which 
would be the Gaussian that, convolved with that of the 
PSF profile, would give a total FWHM = FWHMpsf 
+ 3sTD, where std gives the standard deviation of the 
FWHM distribution measured for all the stars that 
were used to compute the PSF profile. For the closest 
systems we are typically limited to about 10 pc and for 
the furthest systems, 40 pc. 

We could estimate reff using Gaussian fits for about 
88% of the knots in the sample. The profile for the rest 
of the knots could not be fitted by Gaussians, generally 
due to their more extended, diffuse and irregular shape, 
to very steep variations of the local background and/or 
to low S/N of the faintest knots. 

4.2.2. Total size 

The total size of the knots was obtained from the 
F814W images, since this filter is less aff'ected by 
extinction than the other photometric filter in this study 
and therefore is more convenient to measure the radial 
extent of the light of the knots. We derived the edge 
of the knots at the point where the flux, as a function 
of radius, equals the local background flux. In order 
to determine this point, we assumed that the knots 
are circular and measured the surface brightnes s in 
concentric rings, as done in iBastian et al.l (l2005ah . 

The size of the knots has been derived then by 
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Fig. 4. — : Relation between effective radius and total radius 
of the knots for which we could estimate the total size value 
with an uncertainty of less than 50%. Both axes are in log 
scale. The fit for the data is shown at the bottom-right comer. 



fitting the radial profile of the light with a power-law 
of the form Flux oc r"^. This function follows the data 
quite well and coiTesponds to the projected density 
profile of a knot with p oc r'^^^ '^ where p is in units of 



Mq pc (IBastian et al.ll2005al) . For slightly-resolved 



knots the PSF may be dominating the size measured. 
We then estimated the half-light radius of the PSF and 
applied the ratio r^.°' / t^^ to the derived size. 

Due to the complex nature of (U)LIRGs, a large 
fraction of them involved in an interaction process, the 
light profile of the local background of the galaxies is 
usually very steep and irregular Thus, the application 
of this method is somewhat limited. Besides, the size 
of the unresolved knots cannot be derived. We could 
measure reliable sizes, with an uncertainty of less than 
50%, for a third of the detected knots. However, we 
could estimate the total size of the remaining knots (as 
long as they are resolved) by using the relation found 
between r^s and the total size (see Fig.|4]l. In any case, 
given these high uncertainties and that the half-light 
radius is an observable which is more generally used, 
the study of the size of the knots in this paper will be 
focused on reff. The value of the total size of the knots 
will be used for the specific purpose of characterizing 
the mass-radius relation of the knots. 

4.3. Completeness tests 

We have run a series of artificial star tests in order 
to assess the completeness and measurement quality 
of the data set. Some considerations were taken before 
applying the tests: 
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- Since each system has a different local background 
level we had to run a test per image. The same 
happens if we consider different filters. 

- In the central regions of the galaxies the local 
background level is higher and steeper than in 
the outer regions. Hence, different completeness 
levels are expected as the galactocentric distance 
increases. To that end we defined two surface 
brightness regions based on the median value 
across the system: the inner region, corresponding 
to high surface brightness levels, and the outer 
region, corresponding to low surface brightness 
levels. Using this approach, we ran completeness 
tests separately for each region and afterwards 
we weighted each computed value by the number 
of knots within each region. Some works have 
used more varying levels of background (e.g., 7 
inlwhkmore et al. 1999, 3 in Haas et al. 2008), but 
given the angular size of most of the systems in this 
study it is not worth defining such number of levels 
in this study. 

- We have also taken into account the size of the knots, 
since most of them are actually slightly resolved. To 
include this effect we also ran two tests separately. 
First, we added artificial stars with a PSF computed 
using foreground stars. Additionally, we artificially 
broadened that PSF by convolving it with a Gaussian 
that gives a FWHM ~ 1.5 x FWHM of the initial 
PSF. We have then taken the average value of the 
completeness limit in each case. 

The tests were run by adding artificial stars from 
apparent magnitudes of 20 to 30, in bins of 0.2 mag. 
We followed several steps to run a completeness for 
each of the 50 magnitude bins defined. We first added 
1000 artificial stars per bin in a random position inside 
a field that covers the entire system. Poisson noise was 
added and only stars which were more than 10 pixels 
apart were selected to avoid contamination between 
them. Next, we determined in which of the two defined 
regions with different surface brightness (inner and 
outer) each artificial star was located. Finally, artificial 
stars were recovered by performing the photometry in 
the same way as we did for the knots. 

Since the artificial stars were placed in random 
positions, some could reside in the same position as a 
real cluster/star/knot in the original image, increasing 



the measured flux. In order to avoid this, we rejected 
those cases where a recovered magnitude differed from 
the input magnitude by more than a certain value. We 
have taken the value of 0.75 mag because, as seen in 
Figure |5] the systematic uncertainty for faint magni- 
tudes is of the order of 0.4 mag, thus a 2cr rejection 
level seems reasonable. False-positive detections with 
that filtering for the brightest artificial stars are not 
expected, since few real bright objects are encountered 
in the images. The completeness at each magnitude 
level is computed as the number ratio of recovered 
over added artificial stars. All the process was done 
10 times and an average completeness value at each 
magnitude level was computed. An example of the 
result for the images of IRAS 04315-0840 is shown 
in Figure |5] 

The mean difference between input and output 
magnitudes of the artificial stars {Amfsi4w and Amf4i5w 
in Figure ID gives an estimate of the aperture correc- 
tion, which is about -0.3 mag for the F814W filter and 
a value close to -0.25 mag for the F435W filter This is 
consistent with the aperture corrections applied in the 
photometry of the sources. Its standard deviation also 
represents a robust estimate of the actual photometric 
error to be associated with each magnitude bin. 
This error is obviously larger than that computed 
by Poisson statistics only, especially toward the 
fainter magnitudes, since the former consider both the 
systematic and random uncertainties (see gray curve 
in Figure |5]l. Therefore, a more realistic uncertainty 
value for the knots with magnitudes near the 50% 
completeness limit would move up to 0.3-0.4 mag. 

The systems are 90% complete within the intervals 
mfsi4w-'2-3-25 and mf435iy=25-26.5, but due to the 
colors of the regions the F435W completeness limit 
usually shifts to brighter magnitudes. We do not see a 
significant difference in the completeness limits when 
the PSF is convolved for the F435W filter (less than 
0.05 magnitudes), but for the red filter the value of 
Amf^i4w is increased by 0.1 mag when the PSF is 
convolved. Some difference is expected, since as the 
PSF widens we recover less flux by using the same 
aperture. 

4.4. Determination of the luminosity function 

The luminosity function (hereafter LF) quantifies 
the number of objects per luminosity bin. It is usually 
characterized in studies of stars in a cluster, star 
clusters and galaxies. The initial mass function and 
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Fig. 5. — : Result of the completeness test for IRAS 04315-0840. In each case, from left to right, the curve is shown for the inner 
region using the PSF of foreground stars, the outer region with the same PSF, and then the same using the convolved PSF. From top 
to bottom, completeness fraction for filter F814W, input-output of the recovered stars in black and Icr value added and subtracted 
in gray are shown, and the same for filter F435W. The solid vertical line indicates a completeness level of 95%, the dashed line a 
level of 90% and the pointed one a level of 50%. 
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the processes that stars, clusters and galaxies undergo 
once they are formed defines the shape of the LF. 
While the globular cluster LF generally appears to be 
well fitted with a Gaussian function (e.g., IWhitmord 
19971) . for the young knots and clusters observed 
in merging systems and other environments, it is 
well described as a power-law distribution: dN oc 
L^"dL;[ or the equivalent form using magnitudes, 
dN oc lOP^'^dMx, Li and M^ being the luminosity 
and magnitude respectively and N being the number 
of clusters. The slopes in both equations are related 
as a = 2.5/3 + 1. The slope generally lies between 1.8 
and 2. 4 in the optical for systems be t ween 30 and 85 
Mpc (IWhitmore & Schweizer 19951: Sch\yeizer et al. 



1 9961: lMilleretal.l l 1997: Sur ace et al.lll998h 



As we have seen in section 14.21 the knots is this 
study consist generally of cluster complexes, espe- 
cially for galaxies located further than 100 Mpc. Thus, 
the LFs of the knots in this study must not be con- 
fused with those in studies of individual star clusters. 
In section 14.51 we explore how the shape of the LF 
changes due to blending caused by the lack of angular 
resolution to resolve individual entities. 

The shape of the LF has usually been deter- 
mined by fitting an equal sized bin distribution. 
Instead, we used a method described and tested 

based on 



in 



iMaiz Apellaniz & Ubedal (120051) . ___ 

D'Agostino & StephensI (119861) . It is based on using 



bins variable in width, such that every bin contains 
approximately an equal number of knots. Hence, the 
same statistical weight is assign e d to e ach bin and 
biases are minimized. Haas et al.l (l2008h have shown 
that this method is more accurate than those based on 
fitting an equal sized bin distribution. The number of 
bins (Nbi„s) is different for all samples and is related 
to the total number of objects in the samples . Based 



on the prescription in lMaiz Apellaniz & Ubedal (120051) 
we computed the number of bins using the expression: 

2/5 

Nhins - 2A^J + 15, where Nc corresponds to the total 
number of knots. 

The fit was performed on the whole range from the 
90% completeness limit to the brightest knot in the 
sample, thus the slope of the LF might be sensitive 
to an incompleteness of about 10%. Using a higher 
completeness limit cuts off a large number of the knots. 
By correcting for incompleteness at the 90 and 95% 
limits the LF would become steeper (we underestimate 
the slope of the LF if we do not correct for incomplete- 
ness). The typic al value of the un derestimate of the 
slope of the LF in lHaas et aP (l2008h is Ao- < 0. 1 for the 



ACS F435W and F814W filters. We have computed a 
similar value, so we adopt it for this study. 

The uncertainties associated with the fit of the LF 
(shown in Tables and figures) should be considered 
lower limits since they do not take into account other 
sources of error In fact, if the brightest bins are 
dropped in the fit the values obtained for the slope in all 
cases agree to within 0.08 dex. Moving the magnitude 
limit higher by 0.5 mag changes the slope generally by 
up to 0.05 dex. Finally, we have also explored how 
the slope varies when performing a linear regression 
fit to the cumulative luminosi ty distribution ( see the 



third case in the Appendix in iHaas et al.l 120081) . The 



variations here are generally within 0.06 dex. Hence, 
combining all these effects, a more realistic uncer- 
tainty for the slopes of the LF that we have obtained 
would be about 0. 1 dex. 

4.5. Distance dependence of the photometric 
properties 

In the analysis we also consider a selection effect 
that we can encounter when studying the photometric 
properties of the knots in (U)LIRGs located at different 
distances. The classical Malmquist bias is caused by 
the fact that systematically brighter objects are ob- 
served as distance (and volume) increases, as a result 
of a combination of the selection and the intrinsic 
scatter of absolute magnitudes. 

Our sample is affected by this bias, since the sys- 
tems are located at different distances (from 65 to 560 
Mpc). As the system stand further, the faintest knots 
become undetected and we can only see the bright 
ones. Besides, for the same angular resolution the 
linear resolution decreases linearly with distance, and 
therefore the knots that we detect at large distances can 
be associations of knots instead of individual knots as 
those observed in less distant systems. Thus, larger 
sizes are computed. Though this is not a complete 
sample. Figure [T] shows that in general more distant 
galaxies have higher infrared lu minosity. T he same 



tendency is seen if the RBGS (ISanders et al.l l2003h 
are taken. This effect has to be subtracted in order to 
observe reliable intrinsic differences between systems 
at different distances. 

We have divided the sample in three Lir intervals, 
to study the dependence of the photometric properties 
with the total infrared luminosity of the system (see 
section |52]i. To asses how the photometric properties 
depend on the distance of the system a simulation 
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Fig. 6. — : Photometric properties for real (top) and simulated (bottom) data. The median values (between brackets < > throughout 
the paper) of the magnitude, color and r^^ distributions are shown. The slope of the LF for the magnitude distributions is also 
presented (a). Throughout the paper N represents the number of detected knots in any given distribution. 



has been performed. The median distance of the low 
luminosity interval (Ljr < 11.65) corresponds to 99 
Mpc, whereas that of ULIRGs (the high luminosity 
interval) is 258 Mpc. Hence, for this sample ULIRGs 
are located a factor of 2.6 further away than the low 
luminous LIRGs. We have simulated the galaxies 
within the first Lir interval as if they were a factor 
of 2.6 further, obtaining the same median distance 
as in ULIRGs. To that end, the HST images have 
been convolved with a Gaussian at the resolution of 
the PSF at the new distance of each galaxy and re- 
binned in order to have the same "pixel size" at that 
distance. Using foreground stars we have computed 
for each image the factor needed in order to preserve 
the flux, since the flux of the point-like objects must 
be preserved. Then photometry has been carried out as 
in the original images, with circular and in some cases 
polygonal apertures. 

The results of the simulation for the F814W filter 
are presented in Figure |6l The magnitudes and the 
sizes increase by about 1.1 mag (a factor of 2.8 in 
flux) and by a factor of 3, respectively. The number 
of detected objects drops by about a factor of 1.6. 
Therefore, in the simulated image, individual knots are 
associations of knots from the original one. As a result. 



not only do the sizes become larger, but its distribution 
also flattens. The median values of r^ff distribution 
in the simulation are 25 and 78 pc for the real and 
simulated knots, respectively. 

The flattening of the LF (from a - 1.72 to 1.56) 
is also expected, since as the knots are artificially 
grouped due to angular resolution effects the bright 
tail of the LF starts to get more populated. Another 
expected result of the simulation is that the color dis- 
tribution does not change significantly. 

If we consider the spatial distribution of the knots 
(inner and outer knots, with projected distances of less 
and more than 2.5 kpc to the nucleus, respectively), 
some changes are also observed; the ratio of inner to 
outer knots shifts from 0.7 to 0.5. Due to the higher 
surface brightness and steeper background profile in 
the inner regions, when convolving the images to get 
lower linear resolution, a large fraction of inner knots 
are diluted by the local background, in contrast to a 
lower dilution degree in the outermost regions. Be- 
sides, more grouping is expected in the innermost re- 
gions, where there is more crowding, hence also losing 
also a larger number of inner knots when compared 
with those located in the outer regions. 
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Fig. 7. — : Color-magnitude diagrams for the identified 
knots. Top: all knots. The red dashed line indicates the 
median of the distribution and the gray band covers 80 % 
of the total number of knots in each plot, rejecting the 10 % 
at either side of the distribution. No reddening corrections 
have been applied to the magnitudes given. A de-reddening 
vector of Av/=1 mag is shown in the top right corner. Middle: 
the same diagrams for the knots in this study within 2.5 kpc. 
Bottom: same as before but for knots further than 2.5 kpc 
from the closest nucleus. 



5. Results and discussion 

5.1. General properties of the knots 

Owing to the high angular resolution and sensitivity 
of the ACS camera, and to the definition of our large 
sample, covering a wider Lir ra nge than previous 
studies focused only on ULIRGs dSurace et al.lll998 , 
uOOO ). we have detected close to 3000 knots. This 
is more than a factor of ten larger than in previous 
investigations and allows us to carry out a statistical 
study exploring the different physical properties of the 



knots over the entire LIRG and ULIRG luminosity 
range. We have then performed the photometry in both 
F435W and F814W filters for these knots. For this 
study we consider the knots with positive detection 
in both filters. The objects identified as the nuclei 
are not studied here, since many of them are believed 
to be contaminated by an AGN (see table [1] and the 
references regarding the spectral class). 

5.1.1. Magnitudes and colors. Average values 

The knots detected in our sample of luminous in- 
frared galaxies have observed F814W and F435W ab- 
solute magnitudes (uncorrected for internal extinction) 



in the -20 < M, 



FS\4W 



<-9 and -19.5 < M, 



F435W 



<-7 

range, respectively, while colors cover the -1 < 
MF435w-MFii\4w ^ 5 range (Figure |7]l. The corre- 
sponding median magnitudes are < Mfsuw >=-11.96 
and < Mf4-isw >=-10.84. The median color of the 
knots corresponds to < Mf4j,5w-Mfi,i4w > - 1.0, which 
clearly indicates the presence of a young stellar 
population (about half of the knots have colors bluer 
than that value). By contrast, the mean color of 
the diffuse local background light, which traces the 
old stellar population in the host galaxies, is 1.72, 
consistent with the typic al colors measured in spirals 
(B-I^1.80: Euetal.lll993l) . 

A small fraction of the knots (2% of the to- 
tal) show extreme red colors iMf43sw-Mfsi4w =3 
mag or redder), likely tracing regions with very 
high internal obscurati on, usually located in the 



more central regions (lAlonso-Herrero et al.l l2006t 



Garcia-Marin et al.l 120091) . If the stellar population in 
these knots was dominated by young stars (i.e., t~10 
Myr), the extinction could be as high as Av=5-7 mag. 
Heavily obscured young clusters with Ay =5-10 mag 
have also been detected in the optical in less luminous 
interacting systems like the Antennae (e.g., clusters 
Sl_l and2000_l, whose extincti on was computed wi th 
the aid of infrared spectroscopy: Mengel et al.ll2008l) . 



5.1.2. Magnitudes and colors. Radial distribution 

In order to understand the star formation in lumi- 
nous infrared galaxies it is important to investigate 
whether the observed properties of the knots (magni- 
tudes, colors and sizes) do show dependence on the 
galactocentric distance. The spatial distribution of 
the knots is shown in Figure |8] where their sizes are 
related to their absolute magnitude and their color to 
the photometric color Mf4i;,w-Mfg,i4w. About 1/3 of 
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Fig. 8. — : Spatial distribution of the clusters measured on the isophotal map of each system. Each F814W image has been smoothed 
to diminish pixel-to-pixel noise and then avoid spurious contours. The field of view is the same as in Figure[2] The symbol 'X' 
marks the knot identified as the nucleus . From bluer to redder colors, knots are blue, green, orange or red, and the sizes plotted 
depend on their Mf smw (see the legend at the end of the figure). North points up and East to the left. 



19 



IRAS 10173+0828 



IRAS 12112+305 



to 

<1 




10 



IRAS 12116-5615 





10 -10 




IRAS 12540+5708 


20 


■• • /^ 





'■ \(^ 




j v.^ yj 't 


20 




40 





20 

10 

IRAS 13001-2339 



-10 







-^ « S; 


(2^ 


or 


^ 8 


D 


'-''"^^ 


N o 


®r 


/ • 


s \® 




('X 


) r 




X • 


/o/ 


eC 




/ 


. 


\2_-- 


oo 


^^. 




../^ 



40 


Oo 






20 


/r? 


^^^^5r~^ 


• 





:^ 


^/^ 


) 


20 




^^-> ^ 


o 






° #. 





10 5 

IRAS 13428+5608 



-10 




to 

<I 





40 20 


-20 -4 




IRAS 13536+1836 




10 


/'^'^ 


— ^ 





ej^Sfc 


; 


10 







10 20 
IRAS 14348-1447 



-40 



-20 -40 40 20 -20 -40 
IRAS 15206+3342 



< 





10 

IRAS 15250+3609 




-10 


10 




o 


^ 







o 


^ 


^ 




10 






Q 


J 



20 10 

Aa (arcsec) 




Aa (arcsec) 



50 50 

Aa (arcsec) 



Fig. 8. — : - Continued 
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the knots lie within a projected galactocentric radius 
of 2.5 kpc while this fraction increases up to about 2/3 
within a radius of 5 kpc. Thus, a substantial fraction of 
the knots are located close to the nucleus at distances 
of less than 5 kpc. Also, as can be seen in the figure, 
in general they represent the population with redder 
colors, likely indicative of higher internal extinction, 
although an older population can not be excluded. Yet, 
there is also an external population of knots (10.8% 
of the total) at distances of more than 10 kpc and 
up to 40 kpc, located along the tidal tails and at the 
tip of these tails. Most of these external knots are 
blue {Mf4j,sw-Mfsi4w <1), and some (5.7%) are very 
luminous {Mf435w <-12.5), suggesting the presence of 
young massive (M>1O^M0 in young stars) objects in 
the outer parts of these systems. Examples of these 
knots can be seen at the northern tip in IRAS 23128- 
5919 and in IRAS 09022-3675, along the northern tail 
in IRAS 14348-1447, etc. 

Throughout the paper, we identify the inner sample 
of knots as those with projected distances of less than 
2.5 kpc to the closest nucleus, and the outer sample as 
those with distances of at least 2.5 kpc. Both distribu- 
tions have a significantly different magnitude and color 
distribution (see Figure|7|, the inner knots being about 
2 mag brighter in both filters and 0.3 mag redder than 
the outer knots. These two distributions are different, 
and do not come from the same parent distribution 
according to the Kolmogorov-Smirnov test (KS, using 
the idl routine KSTWO; Press et al. 1992). The inner 
knots also have a larger color dispersion with a tail 
toward red colors. This is likely the effect of having 
higher and patchier internal extinction in the inner- 
most regions, in agreeme nt with recent spectroscopic 
studies of these systems lAlonso -Herrero et alj|2006t 
Ibarcia-Marin et al. 2009). Additionally, the brighter 
blue magnitudes in the inner regions suggest that there 
must be more young star formation there than in the 
outer field. 

5.1.3. Effective radius of the knots 



As mentioned in section 14^21 the angular resolution 
hmit (and therefore upper sizes for unresolved knots) 
varies from around 10 to 40 pc, depending on the 
distance to the galaxy. The detected knots do show 
a wide range of sizes ranging from unresolved (12%) 
to a few with very extended sizes (a radius of 200 pc, 
and up to 400 pc). The median size of the resolved 
knots is 32 pc, with about 65% of them smaller than 
40 pc. This is significantly larger than the largest knots 
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Fig. 9. — : Effective radius for inner and outer knots. 
Throughout the rest of the paper the resolved sizes are 
grouped into blank bins. The knots with unresolved sizes 
are grouped into hatched bins, above the blank bins. 



found in other systems, and several times larg er than 



their mean knot radii (e.g., 10 pc in Whitmor e et al. 



1999b . Given the distances of the galaxies in the sam- 
ple (from about 65 Myr to more than 500 Myr) and the 
resolution limit, in general the knots identified in our 
galaxies are likely aggregates of individual clusters, 
which would increase their (apparent) size. According 
to the simulations (see section 4.5), the apparent size 
of the knots changes from an average of 25 pc to 78 pc 
when galaxies at an average distance of 100 Mpc are 
moved at distances 2.6 times further away. 

If we consider the sample of inner and outer knots, 
their distributions are similar (Figure |9ll, so there is 
no evidence of a radial dependence of the sizes at the 
angular resolution of our data. 

5.7.4. The bluest knots. Age and mass estimates 

The determination of ages and masses of stellar 
populations is highly degenerate when using a sin- 
gle color index. In addition, uncertainty (or lack 
of knowledge) in the internal extinction in heavily 
obscured systems like the luminous infrared galax- 
ies considered here, adds more degeneracy. How- 
ever, a first-order estimate of the age and mass can 
be obtained for a specific range in magnitude and 
colors where degeneracy can be minimized. In or- 
der to perform this kind of analysis, two indepen- 
dent synthetic stellar population models are considered 
in this study: the evolutionary synthesis code star- 
burst99 (hereafter SB99) y 5.1 (Leithereret al. 199^; 
Vazquez & Leithereiil2005b and the code bv lMaraston 
(120051). The former is optimized for young population 
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and also models the ionized gas whereas the latter 
gives a rigorous treatment of the thermally pulsing 
asymptotic giant branch (TP-AGB) phase, which is 
relevant within the interval 0.1-1 Gyr In both cases 
we assume instantaneous burst models with a Kroupa 
Initial Mass Function (Kroupa 2002) over the range 
0.1-120 Mq and solar metallicity. These models are 
normalized to 10^ M©. Mass estimates must be mul- 
tiplied by a factor of 1.56 if the Salpeter IMF is 
considered. 

The evolutionary tracks using both SB99 and 
Maraston SEDs are presented in Figure [TOl The shift 
in color at 1 Gyr in the SB99 models is the result 
of using Geneva tracks (optimized for young ages) 
at ages younger than that and Padova tracks at older 
ages, since in completeness sets in for Gen eva models 
at t > 1 Gyr J Vazquez & Leitheren l2005l) . At very 
young ages (log t < 6.5) the nebular continuum, which 
is considered in SB99, reddens the colors by up to 0.6 
mag. From t > 50 Myr (log t = 7.7) up to 1 Gyr the dif- 
ferences between both models also become significant. 
However, regardless of the model used, the degeneracy 
in color within the range Mf43sw-Mffi,i4w =[0.6,1.4] 
is considerable for two reasons: stellar populations 
with ages within the interval [50-500] Myr have the 
same Mf435w-Mfsi4w color in the range 0.6-0.8; the 
Red Super Giant population reddens the tracks at 
ages t~ 7-14 Myr, when this population dominates the 
luminosity, hence making it indistinguishable from 
a population of up to t~l Gyr. Considering this, 
there is a clear age degeneracy for the observed color 
interval Mf4T,5w-Mffii4w =[0.6-1.4] that translates into 
an uncertainty of about a factor of 100 in the mass 
estimates. Hence, although the ignorance of the 
extinction contributes as well, the age degeneracy 
represents the primary source of uncertainty. 

However, as already mentioned in previous sec- 
tions, a non negligible fraction (15%) of all the knots 
in this study are blue, implying young stars and 
some reddening. In particular, knots with colors 
Mf4j,5w-Mfiii4w ^ 0.5 (hereafter blue knots) must 
be young and should not be very much affected by 
degeneracy and extinction, according to the stellar 
population models considered in this study. They are 
bright {<Mf4is\Y >=-l 1.5) and in general found in the 
outer regions of the systems, with some exceptions 
such as IRAS 13536-H1836 and IRAS 20550-H1656 
where they are also detected very close to the nuclei 
(see Figure [8]l. Some first-order estimation of the age 
and mass can be achieved for the blue knots with 
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Fig. 10. — : Color as a function of age (left) and color- 
magnitude track (right) according to the Stellar Synthesis 
Population SB99 models with and without the nebular 
continuum (solid and pointed line, respectively) and to 
the MarastonOS models (dashed line). The curves are 
normalized to 10* Mq. From logt ~ 6.8 on both SB99 curves 
are identical, since the inclusion of the nebular continuum 
is only relevant at younger ages. The de-reddening vector 
(Av=l mag) is shown in both plots. For the whole article it 
is computed using the curve of lCalzetti et al.l ( I200C1) . 

colors Mf4i5w-Mfs\4w ^ 0.5. In that color interval 
the mass uncertainty reduces considerably down to a 
factor of 2-4 if a single stellar population is considered. 
That color can only indicate young population (if 
MF435w-Mff,i4w <0.5, then t< 30 Myr), hence the 
SB99 tracks will be used for this estimation. In case of 
having some degeneracy (contribution from the RSG 
branch), the average age is taken, in order to minimize 
the uncertainty of the estimation. 

Note that the embedded phase of the young popu- 
lation in any environment can last few Myr and that 
(U)LlRGs are known to have young population highly 
enshrouded in dust. Previous works have measured 
very high extinction values for clusters younger than 
about 3 Myr in different environments (Larsen 2010|). 
They have measured extinction values as low as Ay = 
0.5 mag for some clusters older than 3 Myr. Finally, it 
is thought that the embedded phase does not last more 
than 5 Myr, since they have observed clusters older 
than that age which are already extinction-free. Hence, 
the knots for which we are trying to make age and mass 
estimates, though young, are likely to be older than 3 
Myr, probably even older than 5 Myr. 

Under these assumptions and considering all caveats 
above, an estimate of the age and mass of the 
blue knots has been performed, assuming a single 
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Fig. 11. — : Ages and masses estimated for the knots 
with color Mf43jiy-Mf8i4iv < 0.5, where the degeneracy 
in the evolutionary tracks is less severe. The dashed line 
histograms refer to the external knots (projected distance > 
2.5 kpc). 

Stellar population (Figure fTTT i. If small extinction is 
present in these regions, the estimated mass would 
be increased by a small amount (up to a factor of 
1.8 for Ay=l mag). Some blue knots appear to be 
more massive than the most massive old globular 
clusters in the Milky Way, which have masses up to 
10'' Mo (Harris 200 li) . However, most of the blue 
knots have masses in the 10"* to few 10^ M© range, sim- 
ilar to those of giant Hii regions and brightest Young 
Massive Cluste rs in other less luminous interacting 
galaxies (e.g., Bastian et al.1 2005bt Whitmore et al. 



1999l;lMengel et al.l2008l:lKonstantopoulos et al.l2009h 
Complexes of less massive star clusters can also 
constitute these blue knots. 

It is interesting to mention that the blue knots in the 
outermost regions (i.e., d > 2.5 kpc) represent the 80% 
of the total number, while most of the star formation 
is likely to be occuiTing within the central regions of 
the galaxies, as mentioned in section 15.1.21 This sug- 
gests that in general in the outermost regions the dust 
extinction is very low in (U)LIRGs. Hence, the most 
obscured knots must be located close to the nuclei, 
as mentioned in the previous section and in agree- 
ment with recent spectroscopic s t udies of these sys- 
tems JAlonso-Herrero et al.ll2006t iGarcia-Marin et al. 



20091). Part of the difference between the number of 
blue knots in the inner and outer regions may also be 
caused by a stronger disruption process in the central 
regions, as claimed bv lHaas et al.l (l2008h . as well as by 
different histories of star formation of the knots in both 
fields. 
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Fig. 12. — : The mass vs. radius relation for knots in 
(U)LIRGs. The solid red line is a power-law fit to the data, 
with index y. The typical uncertainty associated with the 
data is shown in the right bottom comer. 



The number of blue knots with estimated masses 
higher than 10'' Mq (~ 10%) drops rapidly in the 
external areas (the dashed line histogram in FigurefTTI). 
Hence, the most massive young knots are also located 
close to the nuclei. 

5.7.5. Mass-radius relation of the bluest knots 

We know that in general the knots represent asso- 
ciations of star clusters created in clumps. A single 
giant molecular cloud (GMC) does not usually pro- 
duce a single star cluster, but a complex of star cluster s 
(lElmegreen & Elmegreenll 19831 : Bastian et alJl2005al) . 
Thus, it is interesting to see how these complexes fit 
into the hierarchy of star formation by searching for a 
relation between the mass and the radius of the knots. 
This can be done for the bluest knots, for which an 
estimation of the mass with an uncertainty of less than 
a factor of 4 was achieved (section l5.1.4l) . Their radius 
was derived as explained in section 14.2.21 We could 
obtain reliable measurements of both values for a total 
of 294 knots, which is about a 10% of the total number 
of knots or about 2/3 of the total number of blue knots. 

The surface brightness profile of the bluest knots is 
well fitted by a power-law with index j0=O.37+O.33, 
which coiTesponds to (assuming sphericity) a three 
dimensional density profile of p oc r'^, with jS= 1.37. 
This power-law profile is similar to that observed for 
GMCs assoc iated with high -mass star-forming regions 
(^=1.6±0.3: IPirogovl 120091) . The similarity between 
the profiles in GMCs and the subsample of bluest knots 
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suggests that the amount of luminous material formed 
is proportional to the gas density. Obviously, measure- 
ments of the gas in molecular clouds in (U)LIRGs are 
needed to prove this. 

Although with high dispersion, the bluest knots 
follow a mass-radius relation (see Figure \T2\ . We 
have fitted a function of the form M oc R'', with index 
7=1.91±0.14. The size and the mass of GMCs in the 
Milky Way are also related, M doud ^ ^licud' ^^ ^ ^'~'^' 



sequence of virial equilibrium JSolomon et al.lll987l) 



This relation has also been found fo r extragalactic 



GMCs (e.g., lAshman & Zepfl 1200 It iBastian et al. 
2005al) . Our result is therefore consistent with the 



bluest knots having the same mass-radius relation 
as GMCs in general, contrary to what is measured 



for ind ividual young star clusters (e.g., IBastian et al. 



2005b). An offset between both relations can be used 
to estimate the star formation efficie ncy within the 
cloud, as done in lBastian et al.l (l2005ah . Data from the 
ALMA observatory, which will be able to reach sub- 
arcsecond resolution, would allow this estimation for 
(U)LIRG systems. 

5.1.6. Luminosity Function of the knots in the closest 
systems 

The large number of detected knots allows us to 
study their luminosity function and compare it with 
that for other less luminous and/or interacting galax- 
ies. However, this comparison should be restricted to 
the closest (U)LIRGs, since given the sizes measured 
complexes of star clusters are detected only for galax- 
ies located at distances larger than 100 Mpc with the 
resolution of the ACS. Although many knots detected 
in our closest systems may still consist of complexes 
of star clusters, in order to minimize distance effects 
(see section l431 l we have considered here a subsample 
that comprises the seven systems closer to this distance 
(see specific distances in table[T]i. 

In order to compute the slope of the LF for the 
knots in both the blue and red filters (as explained 
in section I4.4l l we assigned a lower limit cutoff to 
the brightest value of the 90% completeness limit in 
order to ensure the same level of completeness for the 
whole subsample. The single power-law fits to our 
blue and red luminosity functions are not significantly 
different, with slopes of a - 1.95 and a - 1.89 
for the F435W (blue) and the F814W (red) filters, 
respectively (see Figure [T3]l. 

The slopes in this study are in disagreement 




Fig. 13. — : Luminosity functions for the knots identified 
in both the red and the blue filters for the subsample of 
(U)LIRGs closer than 100 Mpc. The line shows the fit up 
to the completeness limit computed for each filter. 



with the much flatter slopes previously measured in 
ULIRGs (Surace et al. 1998). Surace and coworkers 
argued that the flattening of the slope is the artificial 
consequence of not resolving individual clusters, but 
associations of clusters (the bright end is overpop- 
ulated), the same bias that we encounter but more 
severe. In addition, other factors could contribute to 
the flattening of the LF in that study, as well as in other 
works (see table O: (i) the slopes are n ot corrected 



for in completeness (like in NGC 7252 in [Miller et al. 
19971 and NGC 4038/39 in lWhitmore et al.lll999h . (ii) 



the statistics are more limited, since the distribution of 
less than 90 knots is fitted within their completeness 
limit interval, and (iii) ACS imaging is more sensitive 
than WFPC2, and since we are observing closer 
systems, the LF can be measured to intrinsically fainter 
magnitudes. 

The slopes measured for our subsample are how- 
ever similar to those obtained for other less lumi- 
nous interacting galaxies such as in NGC 4038/39, 
NGC 3921, NGC 7252 and Arp 284 (see table E]). 
As the studies start to resolve individual clusters and 
corrections for incompleteness are applied, the slopes 
tend to consolidate around a value of 2-2.2. For 
instance, reft (between 5 or unresolved, and 50 pc) 
derived in early studies dWhitmore & Schweizerll995h 
for clusters in NGC 4038/39, indicate that some blend- 
ing existed and thus the slope of the LF is some- 
what flatter than in more recent studies which detect 
only individual clusters and corrections for incom- 
pleteness are applied as well t Whitmore et alJll999l 
i2010:) . Therefore, the result obtained in this study 
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extends the universality of the slope of the luminosity 
function regardless of the intensity of the star forma- 
tion in interacting galaxies, to extreme star-forming 
systems like the luminous and ultraluminous infrared 
galaxies, at least for systems closer than 100 Mpc. 

5.2. Properties of the knots as a function of Lir 

For most (U)LIRGs the infrared luminosity comes 
mainly from the re-emission by dust of light emitted 
in intense episodes of star formation. In this study we 
span a range of a factor of 15 in infrared luminosi- 
ties from LirH - 11.39 to Lir = 12.54. In order to 
investigate whether the properties of the star-forming 
knots do show dependence on the infrared luminosity 
(i.e., star formation rate), we have divided the initial 
sample into three Lir intervals. We aim at covering the 
low, intermediate and high luminosity regimes while 
preserving a similar number of systems per luminosity 
bin for the subsequent statistical analysis. We have 
therefore chosen the following luminosity intervals: 
Lir < 11.65 (low), 11.65 < Lir < 12.0 (intermediate) 
and Lir > 12.0 (high). The number of systems that lie 
in each interval is 11, 9 and 12, being at an average 
distance of 99, 166 and 258 Mpc, respectively. A 
summary of the properties of the knots per luminosity 
bin is given in Table |4] and will be discussed in detail 
in subsequent subsections. 

As shown in Figure[Tl systems with higher infrared 
luminosity are in general intrinsically more distant. 
We have tried to get rid of any distance effects by 
selecting a subsample of systems located at a similar 
distance (those at 65-75 Mpc) and sampling all the 
Lir bins. This subsample, not homogeneously dis- 
tributed, includes only one system in the intermediate 
and high luminosity bins (4,1 and 1 systems, respec- 
tively per luminosity bin). Besides, the results based 
only on the ULIRG Arp 220 (IRAS 15327+2340) are 
not reliable since it is known to suffer extreme ob - 
scuration from its silicate features JSpoon et al.ll2007l) . 
In order to achieve more reliable statistics we need 
to include a larger number of systems in the defined 
luminosity bins and hence we inevitably have to deal 
with distance effects. 

Finally, for each Lir bin, we have excluded the 
knots of systems located at the two extremes of the dis- 
tance scale, in order to diminish distance effects within 
the intervals. Thus, the statistics that are shown first in 



^We remind the reader that for simplicity, we identify the infrared 
luminosity as Lir (Lq ) = log (L/r[8 - lOOOpm]) . 



Table |4]for the intermediate luminosity interval do not 
take into account IRAS 04315-0840, and for the high 
luminosity interval the knots in IRAS 15206-1-3342 and 
in IRAS 15327-1-2340 are not considered either. How- 
ever, for completeness, the statistics is also provided 
considering all the sources in each Lir bin. 

5.2.7. Magnitudes and colors of the knots as a func- 
tion ofLiji 

Knots show a dramatic increase (x 12) in their blue 
and red luminosities (i.e., 2.7 mag) with the infrared 
luminosity of the entire galaxy while basically pre- 
serving their colors (see histograms in Figure [HI and 
average properties listed in Table |4|. Our simulations 
(see Figure |6] and section 14.5b suggest that distance 
effects can explain only part of the differences revealed 
by the data. If the same type of galaxies are located 
at the average distance of the low, intermediate and 
high luminosity subsamples, the distance effect does 
not change significantly the colors while increasing 
the luminosity (absolute magnitude) by only a factor 
of 2.8 (1.1 mag). Therefore, we conclude that knots 
in highly luminous galaxies are at least a factor of 
four more luminous than knots located in systems with 
lower infrared luminosity. 

It is interesting to mention that when compar- 
ing these results with those independently obtained 
for nearby, less luminous systems (the weakly and 
strongly interacting galaxies M5 1 and NGC 4038/4039, 
respectively), and more distant ULIRGs (Surace et al. 
1998), the same trend with infrared luminosity appears 
(see Figure [15] and table |4]l. In nearby systems like 
M51 and NGC 4038/4039, the angular resolution is 
such that individual clusters, instead of aggregates as 
in more distant galaxies, are detected, and therefore 
the difference in the luminosity of clusters and knots 
can be partly due to a distance effect. Knots in distant 
ULIRGs tend to be as luminous as knots detected in 
our closer high-luminosity galaxies (also ULIRGs) 
while having colors similar to those knots detected in 
our sample and clusters in M51 and NGC 4038/4039 
(The Antennae). The most luminous clusters in 
the Antennae cover a luminosity range close to the 
median value for knots detected in our low-luminosity 
systems. Finally, clusters in the very low luminosity, 
weakly interacting M5 1 galaxy are on average about 3 
magnitudes fainter (i.e., xl5 less luminous) than knots 
in our low-luminosity galaxies. 

If we assume a statistically similar extinction and 
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Table 3 
Slopes of the LF computed for this work and for other interacting systems 



Galaxy 


Ol " (Mpc) 


Band 


Slope LF (af 


Uncertainty 


Ref. 


ULIRGs Surace 


422 


I 


1.39 


0.08 


[1] 


NGC 3921 


84.5 


V 


2.10 


0.30 


[2] 


Nearby (U)LIRGs 


72.8 


B,I 


1.95,1.89 


0.02,0.04 


This work* 


NCG 7252 


62.2 


V 


1.84 


0.06 


[3] 


Arp 284 


33.6 


I 


2.30 


0.30 


[4] 


NGC 4038/39 


27.5 


V 


1.78 


0.05 


[5] 


NGC 4038/39'-' 


27.5 


V 


2.01 


0.11 


[6] 


NGC 4038/39 


27.5 


V 


2.13 


0.07 


[7]* 



^In some works they also fit a double power-law, but the slopes given in this table correspond 
to a single power-law fitting. 

''D istances taken from NED. For the (U)LIRG samples (this work and that bv lSurace et alj 
1 19981) the mean distance is given. 

'Computed for cluster-rich regions on the PC chip. 



1997h . 
I99I i 



Referenc es. — [11 ISurace et all h998l), [21 ISchweizer et"ai] jl996h , [31 iMiller et al 

[41 IPeterson et al.' ('2009^. [51 'Whit more & Schweizerl h995h . [61 IWhitmore et~ 

and [7 1 Whitmore et al. (2010). In asterisk, studies that have used the prescription 



in lMai'z ApelTaniz & Ubedal J2005h . 



Table 4 
Photometric properties of the sample as a function of Lir 



System/s 
(1) 



(or intei'val) 

(2) 



Number of 
systems 

(3) 



knots per 

system 

(4) 



<Di> 

(Mpc) 
(5) 



< Mnuw > < Mpiiiw > <C> < r'J 



(6) 



(7) 



(8) 



(pc) 
(9) 



o-LF 

F814W 

(10) 



ffLF 

F435W 
(11) 



Lir < 11.65 (all) 
1.65 <L,R < 12.{ 



(U)LlRGs 
(this work) 



.^-^ J MR < 12.0 
Lir > 12.0 
1 1.65 < Lir < 12.0 (all) 
LiR> 12.0 (all) 



11 


113 ± 102 


99 ±26 


-11.36 


-10.13 


1. 13 


25 


1.80 ±0.02 


1.77 ±0.02 


8 


84 ±59 


178 ±27 


-12.74 


-11.84 


0.81 


55 


1.73 ±0.04 


1.79 ±0.04 


10 


38 ± 15 


246 ± 72 


-14.09 


-12.97 


0.97 


81 


1.83 ±0.03 


1.78 ±0.04 


9 


115 ± 109 


166 ±45 


-12.11 


-11.18 


0.82 


45 


1.78 ±0.04 


1.78 ±0.04 


12 


47 ±41 


258 ± 125 


-13.39 


-12.20 


1. 10 


59 


1.78 ±0.04 


1.78 ±0.04 



M51 

NGC 4038/4039 

ULIRGs (Surace98) 



10 

10.7 
Lir > 12.0 



881 10.6 

~ 10* 27.5 

12 ±10 520 ±210 



-8.4 
-11.7 
-14.7 



-7.6 
-11.0 

-13.7 



0.7 
0.9 
0.8 



16.8 
65 



Note. — The first three rows show the statistics for the systems considered in the text (some few have been excluded in the intermediate and high infrared 
luminosity intervals; see beginning of Section l5.2l The following two rows implement the statis tics for a ll the systems that define the Ljr interval. Below 
the horizontal middle line the v alues for clusters and knots from other works are shown (M51. lBik et a l. 2003; NGC 4038/4039, Whitmore et al. 1999J; 
and ULIRGs. ISurace et al.iri998l) . Col (1): galaxy name or designation for a sample of galaxies. Col (2): infrared luminosity (or interval) of the 
galaxy/sample. Col (3): number of systems that comprises each sample. Col (4): average number of knots per system. Note that even though more 
than 10 clusters were detected, the values for NGC 4038/4039 refer to the 86 brightest clusters. Col (5): luminosity distance or average value. Col (6): 
median value of the Mfgi4w absolute magnitude distribution. Col (7): median value of the Mf^^^^y absolute magnitude distribution. Col (8): median 
value of the Mf4jj\Y-Mfg[4w color distribution. Col (9): median value of the etfective radius distribution. Col (10): slope of the F814W luminosity 
function distribution. The slopes of the LFs for interacting systems in other works are aheady showed in Table [3] Col (11): slope of the F435W 
luminosity function distribution. 
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Fig. 14. — : Photometric properties of the sample in three Ljr intervals, excluding 3 systems in the sample, as explained in 
the text. The dashed vertical line indicates the median value of the distributions of magnitudes Mf43sw and Mf smw , and color 

Mf435w-Mfg[4W 



age for the star-forming knots regardless of the Lir of 
the system, the most plausible explanation for the 
luminosity excess measured in the intermediate and 
high luminosity systems has to invoke a mass and/or 
density effects: 

• Mass effect. More high-mass knots are sampled 
in the most luminous systems. 

• Density effect. The knots are aggregates of an 
intrinsically larger number of clusters as the 
infrared luminosity of the system increases. 

The star formation rate as well as the gas con- 
tent in ULIRGs is higher than in less luminous sys- 
tems and as a consequence, these systems form more 
clusters. Therefore, it is matter of simple statistics 
rather than a difference in physical formation mech- 
anisms that in ULIRGs we find the brightest knots. 
Given the stability of the colors among the different 
infrared luminosity sample, they likely represent the 
most massive knots. Thi s stochastic process, k nown 
as size-of-sample effect (IWhitmore et al.l l2007h . can 
explain the demographics of star cluster systems in 
the merging environment, where the star formation is 



enhanced. Therefore, although density effects cannot 
be discarded our result is compatible with a mass effect 
(which in turn corresponds to a size-of-sample effect), 
explained by statistics. 

5.2.2. Spatial distribution of the knots as a function 
ofLi^ 

One important aspect to investigate is whether the 
spatial distribution of the knots shows dependence on 
the luminosity of the system. We measure a ratio 
of inner to outer knots i^hmerl^ outer) that changes 
with Lir from 0.68 to 0.49 and 0.31 for the three 
(low-, intermediate- and high-) luminosity intervals, 
respectively). According to our simulations, these 
differences are consistent with being purely due to a 
distance effect (see table |4] for the average distance 
of the different luminosity bins), and therefore no 
evidence for the dependence of the spatial distribution 
of the knots on the infrared luminosity of the systems 
is found. 
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Fig. 15. — : Color-magnitude diagrams of all the identified knots (black) at different Lir intervals, compared with other interacting 
systems from the literature. The gray band covers the 80 % of the total nu mber of knots in each plot, as in Figure[7] The colors refer 
to clusters/knots measured in M51 (orange) taken from lBik et al.l ( 120030 ■ in the ULIRG sample of ISurace et al.l 1 11998 ) (pink) and 
in NGC 40 38/4039 (blue, those y ounger than 10 Myr; green, those between 250 and 1000 Myr; and red, those older than 1 Gyr), 
taken from lWhitmore et alj ( Il999l) . 



5.2.3. Effective radius of the knots as a function of 

If there is a change in the luminosity of the knots 
with the bolometric luminosity of the system, their 
size could show a variation as well. At first look 
this seems to be the case (see Figure [T4l and [TSIl. As 
the luminosity of the system increases, the distribu- 
tion of Teff shows a broadening with a non negligible 
fraction of the knots having sizes larger than 100 pc. 
This is more clear for the high-luminosity systems 
(ULIRGs). However, as our simulations have shown 
(see Figure|6]l, this changes can be explained purely as 
a distance effect. On the one hand, reff distribution of 
knots in ULIRGs is similar to the simulated r^s distri- 
bution. Moreover, the median value of rgff in ULIRGs 
is 81 pc (see table |4|, slightly larger, but consistent 
within the uncertainties, with the 75 pc derived from 
the simulations. In fact, the sizes of the knots in the 
closest systems (i.e., Dl <100 Mpc), where individual 
clusters are almost resolved, are closer to the sizes of 
clusters measured in less luminous interacting galaxies 



(e.g., see clusters of NGC 4038/4039 and NGC 7252 
in Fig [T6b . Therefore, with the present angular reso- 
lution there is no evidence for a variation in size with 
luminosity. 

5.2.4. Mass-radius relation of the bluest knots as a 
function of L\s^ 

We have seen that in general terms the bluest knots 
in our sample follow a mass-radius relation with index 
y ~2, similar to that in GMCs and complexes of 
clusters in less luminous interacting systems (see sec- 
tion |5.1.5l l. Here we have the opportunity to character- 
ize this relation as a function of the infrared luminosity 
of the system. 

Though with high dispersion, the coiTelation holds 
for the low and intermediate infrared luminosity inter- 
vals (see Figure [T7] i. The fit for the bluest knots in 
ULIRGs gives a shallower index, though it could still 
be consistent with a value of y-1.9 given the large 
dispersion on the index values. The fit is affected 
in any case by the small number of blue knots in 
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NGC 7252 
ULIRG Surace98 







1000 



Fig. 16. — : Absolute F814W magnitude (uncorrected for internal extinction) as a function of size for clusters in interacting systems. 
Note that the x-axis is in log scale. We have divided the knots in our sample in three groups dep ending on their infrared luminosity. 
Posit ive values for clusters in the Antennae were taken from Whitmore & Sc hweizen ( Il995l) . in NGC 7252 (with Ljr =10.70) 
from lWhitmore et alj (Il993i) and in the ULIRG sample from iSurace et al. til998i1 . The inset plot shows the same for the unresolved 
knots in this study. 



ULIRGs and probably by distance effects. Therefore, 
according to our data, there is no clear variation of the 
mass-radius relation with the infrared luminosity of the 
system. 

5.2.5. Luminosity Function of the knots as a function 
ofLiR 

The derived slopes of the luminosity functions with 
values between 1.7 and 1.8 (see tabled do not show 
any clear trend of variation with the luminosity of 
the system. However, our simulations of the dis- 
tance effect show that the slope of the LF becomes 
shallower by about 0.15 dex when moving the low- 
luminosity systems to the average distance of the high- 
luminosity subsample (see Figure |6]l. This change in 
the slope is due to the reduction in the linear resolution 
such that knots that appear resolved in the closer, 
low-luminosity systems appear to be as more lumi- 
nous aggregates in the more distant luminous systems. 
Therefore the distance-corrected LF slopes for the 
intermediate and high-luminosity subsamples would 
be about 1 .9-2, in the range of those measured in other 
low-z interacting galaxies (see tableO. 



5.3. Properties of the knots as a function of inter- 
action phase 

Strong interactions and mergers provide a natu- 
ral laboratory for probing how star formation is af- 
fected by tidal forces and major rearrangements of the 
gaseous component during the different phases of the 
interaction. Most of the systems in our sample are 
interacting systems, hence we have divided the sample 
into four groups, depending on the morphology seen in 
the F814W HST images. Each group represents snap- 
shots of the different phases of the interaction/merger: 
first approach, pre-merger, merger and post-merger 
The morphology class is explained in section 12.21 and 
the selection is shown in Table[T] In order to know how 
the star formation evolves with the interaction process, 
in this section we study the photometric properties of 
the knots for each of the four morphology groups. 

5.3.1. Magnitudes and colors as a function of inter- 
action phase 

The magnitudes and colors of the knots in the first 
three interaction phases differ significantly from those 
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11.65 <L|g< 12 

— Y=2.12±0.23 



H 1 h- 




or an older population. In the following we explore 
both alternatives. 



1.5 2.0 2.5 

log (radius/pc) 

Fig. 17. — : The mass vs. radius relation for laiots in 
(U)LIRGs as a function of Ljr. The solid red line is a power- 
law fit to the data, with index y. 



of the knots in the post-merger phase (see results 
in Table |5]l. Even if in a strictly statistical sense 
(KolmogorofF-Smirnov test at a confidence level of at 
least 98%) the population of knots in the first three 
phases are different, their median values and distri- 
butions are very similar suggesting the lack of strong 
differences. This is true in particular if, in order to 
minimize the distance effect, only galaxies at similar 
distances in each phase are considered (see table |5] 
and Figure [TSll. However, the population of knots in 
the post-merger phase are significantly more luminous 
(median I-band magnitude difference of up to 2 mag), 
and significantly redder than in any of the three ear- 
lier phases. Besides, the population of knots in the 
post-merger phase shows in its color distribution an 
extension toward redder colors (> 3 mag). The fact 
that knots in the post-merger phase have redder colors 
means that they likely represent either a more obscured 



1. 



In order to investigate whether the knots in 
the post-merger phase are intrinsically more 
obscured than in any of the previous phases, 
we define the infrared to ultraviolet luminosity 
ratio (IR/UV) as 



IRjUV = 



^IR 



AL^FUV) + AL^iNUV) 



where Lir stands for the infrared luminosity, A 
to the pivot wavelength of a given filter and L^ 
to the flux in the n ear (NUV) and far (FUV) 
ultraviolet emission (IHowell et alJl2010h . This 
ratio is a useful empirical measurement of the 
dust absorption and emission in these heavily 
obscured systems. We searched for UV images 
in the Galaxy Evolution Explorer (GALEX) cat- 
alog, using the Multimission Archive at STScI 
(MASlQ) facility. With a typical angular resolu- 
tion of about 4-5" in the FUV-NUV (A = 1528- 
2271 A) range, we can estimate an average value 
of IR/UV for each system assuming that the far 
infrared and ultraviolet flux is emitted in the 
same regions. Even if this is not necessarily true, 
taking into account that most of the nuclear star 
formation is obscured in (U)LIRGs, we can have 
a first-order estimate on which systems have a 
higher average extinction. 

A total of 16 systems were detected in the GR5 
data set, most of them in the All-Sky Imaging 
survey (ASO. The other systems were either 
not detected or not observed. In this dataset, the 
images are akeady reduced, the sources detected 
and the photometry performed. In general the 
photometry was performed using an aperture 
that covers the entire galaxy (if two galaxies 
form the system we add up the flux values). We 
then take the value of the flux provided by the 
data set facility. 

With the flux in the far and near UV and the 
infrared luminosity we have estimated an aver- 
age value of the IR/UV ratio for each system. 
Although the assumption by Kennicutt (19981) 



''http://archive.stsci.edu/ 

'' Visit the GALEX website, [iittp://w 



information 



v.galex.caltech.edii/] for more 
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Fig. 18. — : Photometric properties as a function of merging phase for systems with similar D^. Dashed vertical lines indicate the 
median value of the distributions of magnitudes M/j8i4w and M/r435iv , and color Mf43jiy-Mf8i4iv , as in Figure [14] The relevant 
parameters of the plots are shown in the first four rows in Table |5] 
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Table 5 
Photometric properties of the sample as a function of interaction phase 



Interaction phase 


Number of 
systems 


knots per 
system 


<Dl > 

(Mpc) 


< ^FSl-iW > 


< Mf4351v > 


<C> 


(pc) 


a LF F814W 


a LF F435W 


First Approach 


3 


153 ± 127 


121 ±49 


-11.25 


-10.26 


0.93 


24 


1.50 ±0.09 


1.63 ±0.06 


Pre-merger 


8 


145 ±81 


126 ±42 


-11.91 


-10.77 


1.06 


35 


1.84 ±0.03 


1.86 ±0.04 


Merger 


6 


139 ±123 


107 ± 37 


-11.23 


-10.24 


0.96 


25 


2.00 ± 0.04 


1.86 ±0.03 


Post-Merger 


5 


34 ±19 


145 ± 50 


-13.35 


-11.49 


1.49 


44 


1.76 ±0.12 


1.58 ±0.17 


First Approach 


4 


131 ±113 


140 ± 55 


-11.54 


-10.66 


0.92 


27 


1.53 ±0.10 


1.67 ±0.02 


Pre-merger 


13 


103 ± 83 


192 ± 98 


-12.29 


-11.12 


1.04 


40 


1.93 ±0.03 


1.86 ±0.03 


Merger 


8 


112±115 


142 ±71 


-11.41 


-10.38 


0.96 


26 


1.83 ±0.02 


1.78 ± 0.03 


Post-Merger 


7 


29 ±18 


213 ±159 


-13.51 


-11.73 


1.55 


49 







Note. — The interaction groups above the horizontal middle hne includes only systems having similar distance. Below the horizontal middle 
line all the systems have been considered, regardless their distance (note the higher dispersion in the distance of the systems here). Apart 
from the first column, which refers to the different interaction phases defined for this study, the rest indicate the same as shown in Table|4] 



Table 6 
Infrared to ultraviolet luminosity ratio as a function of interaction phase 

Interaction phase N,,,, < IRAJV > IR/UV < IR25AJV > IRjs AJV 

range range 

(1) (2) (3) (4) (5) (6) 



First Approach 


2 


- 


[15,130] 


- 


[4,11] 


Pre-merger 


6 


48 


[14,189] 


7 


[2,45] 


Merger 


6 


248 


[114,1235] 


32 


[14,117] 


Post-merger 


2 


- 


[385,895] 


- 


[86,93] 



Note. — Since the IR/UV ratio is not affected by distance effects (we compute it for 
the whole system) we give the statistics for all the systems per interaction phase with 
available UV image. Col (1): interaction phase. Col (2): number of systems with IR 
and UV data available. Col (3): median of the IR/UV ratio when more than two values 
are derived. Col (4): range covered of the ratio IR/UV for each interaction phase. Col 
(5): same as (3), but in this case only the contribution of the flux at 25pm is taken to 
compute LiR. Col (6): same as (4), but only the contribution of the flux at 25yum is 
taken to compute Lir. 
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that the great majority of the luminosity from 
young stars would be absorbed by dust and re- 
radiated in the far-IR is likely to be correct, ob- 
servationally there may be other contributions to 
the total IR luminosity from older populations of 
stars or other luminosity sources (e.g., see the re- 
view by lTuffs & Popescull2005h . For that reason 
we have also computed the IR/UV ratio taking 
only the contribution at 25/vm in IRAS (IR25), 
based on the good relation between the mid- 
IR 24jum-25yum (MIPS -IRAS) and the number 
of ionizing photons as derived from extinction 



corrected Pa^ and Hg lurninosi t ies (IWu et al. 



2005 



Alonso-Herrero et alJl2006l: ICalzetti et al. 



20071). This way we select only hot dust, heated 
by the radiation of young population, and avoid 
any contribution from colder dust, heated by 
older stars. 

Bearing in mind that we are dealing with a small 
number of statistics, there is a trend toward 
higher values of the IR/UV ratio in not only 
the post-merger but also the merger phase, with 
respect to early phases represented by the first 
approach and pre-merger (see values in Table|6l). 
Although the absolute values are different and 
the refinement is better when using IR25, the 
trend is practically the same if we use one in- 
frared luminosity or the other. This is not sur- 
prising, as IRAS measurements poorly sample 
the output of the cold dust (the longest band is 
at lOOyum). Since the IR/UV ratio is related to 
an average dust absorption, our result suggests 
that knots in the late phases of an interaction 
could be suffering on average intrinsically more 
extinction than in early phases. For instance, 
an additional average extinction of Ay =1 can 
already redden the colors by about 0.6, accord- 
ing to the SB99 models used in section |5. 1 .41 
Were this the case, the population of clusters 
in the late phases of a merger would be more 
luminous than in previous phases. Therefore, 
we would need to invoke a physical process 
responsible for this significant increase in lumi- 
nosity. Note however that dust absorption shows 
signific ant structure on scales o f a kiloparsec 
or less dGarcia-Marin et al.ll2009h . and therefore 
the measured IR/UV ratios do not necessarily 
represent the intrinsic absorption toward the in- 
dividual optically-selected clusters. A larger 
set of high angular resolution UV and mid-IR 



imaging is required before any firm conclusion 
is derived. 

2. An alternative scenario to explain the differ- 
ences in color and luminosity of the knots iden- 
tified in the post-merger phase with respect to 
early phases c an invok e evolu tionary effects. As 
proposed by iKroupal ( 119981) . many dozens or 
even hundreds of massive individual star clus- 
ters with masses of lO^-lO^M© could merge into 
extremely massive super star clusters (SSCs) 
after a few hundred Myr This hierarchical 
star formation evolution has also been suggested 
to be happening in complexes of clusters in 
M51 JBastian et al.l E005a). As the interaction 
proceeds, the smaller SSCs would be destroyed 
by disruptio n effects (e.g., by the mechanisms 
proposed in IWhitmore et al.l l2007b and a few 
massive individual clusters would still be able 
to form. Note that, although with a large dis- 
persion, while the number of knots per system 
stays about the same (140-150) in the first three 
phases, this number drops by a factor of 4 to 
5 in the post-merger phase (34 knots per sys- 
tem). Thus, in the post-merger phase systems 
we might be detecting only the evolved, ex- 
tremely massive SSCs and a few massive in- 
dividual clusters. The knots with the reddest 
colors could be evolved SSCs as massive as 
10^"'' Mo with an optical extinction (Ay) of 
3 mag, or massive young clusters with higher 
extinction. 



Kroupal (Il998h estimated that after 100 Myr 
or so an object with these characteristics can 
be formed. If the merging of clusters started 
during the early phases of the interaction we 
would have already detected some of them in 
the merger phase. Given the degeneracy of 
the Mf435v(?-Mf 8i45y color we could only set an 
upper-limit to the number of merged superclus- 
ters in this phase. Thus, with a color higher than 
Mf43$w-Mfsi4w - 0.6 and a luminosity higher 
thanMf8i4vi/ = -15 an object could be more mas- 
sive than 10^ Mq, should it be around or older 
than 100 Myr. About three such knots per sys- 
tem are detected in our sample of (U)LIRGs. 

Obviously, photometry at other wavelengths 
is necessary to better constrain the ages and 
masses of these knots and thus to check the 
vahdity of this scenario for (U)LIRG systems. 
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Likewise, simulations of galactic encounters 
help us know if this scenario can occur in 
this environment. A detailed comparison of 
the present observational results with state- 
of-the -art high-resolution n umerical simula- 
tions (IBournaud et al.ll2008bl) will be presented 
elsewhere (Miralles-Caballero et al. 2011, in 
preparation) to further investigate this scenario. 

5.3.2. Spatial distribution as a function of interac- 
tion phase 

In order to know if major mergers have an impact 
on the spatial distribution of the knots, we have in- 
vestigated the relative fraction of inner (inside 2.5 kpc 
radius) and outer (outside 2.5 kpc radius) knots in the 
four interaction phases defined above. Although the 
number of systems in each phase is small (3, 8, 6 and 
5 for the first approach, pre-merger, merger, and post- 
merger) the ratio Ni„„er/Noui<T is 0.50,0.65,0.52,0.68 
from the first approach to the post-merger. 

We do not see any correlation toward a more dis- 
persed or concentrated spatial distribution of the knots 
with interaction phase. However, in the pre- and 
post-merger phases, where the magnitudes are higher 
compared to the other phases, there is some evidence 
of greater concentration toward the nuclear regions. 
The greater the number of knots in the central re- 
gions, the harder it is to detect them, and therefore, 
only the brightest ones are identified given the local 
background conditions (bright and steep light profiles). 
Hence, this effect can contribute to a small increase in 
the luminosity distribution of the knots in these phases. 
Assuming that the knots in the pre-merger phase have 
a luminosity similar to that in the first contact and 
merger phases, this effect would increase the median 
luminosity by a factor of almost 2. Since the ratio 
^iimer/^outer in this phasc is similar to the ratio in the 
post-merger phase, an additional increase of a factor of 
about 5 would still be unexplained by the first scenario 
explored in the previous section. A larger sample of 
galaxies would be required to establish a firmer result. 

5.3.3. Effective radius as a function of interaction 
phase 

All size distributions are statistically different with 
a high degree of confidence (higher than 99%), yet 
the median value of r^ff is similar during the three 
first interaction phases indicating that these differences 
must be rather small. However, knots in the post- 
merger phase tend to be larger by factors of 1.3 to 2 



than in any of the previous phases (see table |5]l. This 
result is not very robust, given the uncertainties in 
the determination of sizes, but if confirmed, it would 
support t he hiera r chical star formation scenario pro- 
posed bv iKroupal (11998 !). According to this scenario, 
the merged super clusters would have a half-mass 
radius of about r/,=45-95 pc. Under the assumption 
of virial equilibrium and isotropic velocity distribution 
the effective and half-mass radius can be related as 
'"eft =3/4r/,. Then, the merged extremely massive super 
clusters would have r^ff =30-70 pc, similar to values 
measured for knots in the post-merger sample. 

5.3.4. Mass-radius relation of the bluest knots as a 
function of interaction phase 

The broad-band emission of the knots for which 
we have measured the mass is very much likely to 
come from young population (i.e., t<30 Myr; see 
section 15.1.41 ). Based on equivalent width measure- 
ments, iBastiatLSLalJ (l2Q05a|) determined the age of the 
complexes in M51 to be in the range 5-8 Myr. On the 
other hand, the phases of interaction are separated by 
hundreds of Myr. Thus, although we cannot study how 
the mass-radius relation of the knots behaves as they 
evolve, we can see if they all have a similar relation at 
a short time after they are born during the interaction 
process. 

The mass-size relation of the bluest knots in the 
three interaction phases for which a fit could be 
achieved is consistent with an index of y=2 (Fig.[T9b. 
Although the line-fit seems to become shallower as the 
interaction proceeds from the first approach up to the 
merger phase (from 7=2.11 to 1.71), the dispersion is 
too high to accept that it varies. Actually, the fit that 
we obtained using all the knots gives a slope inbetween 
both values (1.91; see section is". 1.5b . 

5.3.5. Luminosity function as a function of interac- 
tion phase 

There appears to be a trend in the slopes of the 
LFs of the knots with interaction phase such that they 
become steeper as the interaction evolves from first 
contact to pre-merger and up to merger and post- 
merger (see table |5]). This is particularly more visible 
when only systems at similar distances and /-band 
LFs, less affected by extinction effects, are considered 
(slopes of 1 .50, 1 .84, 2.0 and 1 .76, respectively). Note, 
however, that the slope of the LF for the post-merger 
phase is computed with less than 100 knots (brighter 
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Fig. 19. — : The mass vs. radius relation for knots in (U)LlRGs as a function of the interaction phase. The solid red line is a 
power-law fit to the data, with index y. 



than the completeness limit), hence it is statistically 
less significant than the other slopes. Also note that the 
systematic uncertainties of the slope of the LF can be 
as high as about 0. 1 dex (the tables give the uncertainty 
related only to the fit), as mentioned in section 14.41 
Thus, the observed trend is very weak. 

If we select only systems closer than 100 Mpc (see 
section I5.1.6I I, the phases are also sampled. We can 
try to further minimize distance effects by computing 
the LF only for these systems. While 400-600 mea- 
surements can be used in each of the other phases, on 
the post-merger phase we are instead restricted to only 
about 30 knots and thus, we will not include this phase 
in this analysis. The slopes derived for the closest 
systems are 1.68-1.68-1.92 (1.75-1.77-1.99) for the 
/- (B-) band LF for the first contact, pre-merger and 
merger phases, respectively. It is interesting to note 
that in this case the LF still becomes steeper when 
the systems evolve from the pre-merger to the merger 
phases. 

(U)LIRGs are known to have high SFRs at all 
phases of interaction. However, this can be con- 
centrated at small scales or more dispersed and, at 
the same time, less or more hidden by dust. For 



instance, 80% of the total infrared luminosity in IRAS 
20550-1-1656 comes from an extremely compact, red 
source not associated with the nuclei o f the merg- 
ing g alaxies but with a buried starburst (Inami et al. 



20101) . If, on average, we detect more young knots 



in the optical at early phases of interaction than in 
more evolved phases, they can overpopulate the bright 
end of the luminosity function, and thus the mea- 
sured slope would be shallower. A forthcoming study 
(Miralles-Caballero et al. 2011, in preparation), com- 
paring the observed L Fs with those computed from 
numerical simulations (IBournaud et al.l |2p08b) (i.e., 
mass function) and stellar population synthesis models 
(i.e., mass to light ratios), will help us learn if the 
observed trend (change of the slope as the interaction 
evolves) is real. Should it be confirmed, then we will 
be able to investigate its origin. 

5.4. Star-forming knots in (U)LIRGs versus star- 
forming clumps in high-z galaxies 

(U)LIRGs represent the extreme cases of low- 
redshift star-forming galaxies associated with spirals, 
interacting galaxies and mergers. Deep cosmolog- 
ical imaging and spectroscopic surveys at redshifts 
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Table 7 
Comparison of knots in ULIGRs with clumps in the high-z Universe. 



System /field 


rest-frame 


redshift 


N clumps 


Mass range 


Total size 


Ref. 




band /line 




per galaxy 


(Mo) 


(kpc) 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


UDF field 


B 


0.5-4 


~5 


10^ - 10' 


1 -2.5 


[1], [2] 


Lensed galaxies 


H„ 


1.7-3.1 


- 


6xlO'*-3xlO'' 


0.3-1 


[3] 


SF galaxies at z~2 


optical g, UV, H„ 


2.2 - 2.4 


4.5 


9x10' -9x10' 


0.3-1.5 


[4] 


GEMS & GOODS fields 


B 


0.1- 1.4 


~4 


5x10"- 10** 


0.2-7 


[5] 


Simulated ULlRGs at z=l 


1 


1 


6.7 


2x10^-7x10' 


0.4 - 3.5 


This work 


TDG candidates in (U)LIRGs 


I 


0-0.1 


~ 1 


10«-4xl0'' 


1 -2 


[6] 



Note. — Col (1): type of systems and deep fields. Col (2): rest-frame photometric band. Col (3): redshift interval of 
the galaxies. Col (4): average number of clumps per galaxy. Col (5): total stellar mass range of the clumps, except for 
the lensed galaxies, where the dynamical mass of the clumps is provided. Col (6): total diameter of the clumps. For 
star-forming galaxies at z = 2 the full-width at half-maximum from the clump radial li ght profiles after su btra ction of the 
local b ack ground is given. Col (7) : refere nces- [1] Elmegreen & Elmegreen (20051); [21 lElmegreerj <2007r) : [31 Ijones et al.l 
bold) : [4l iF6rster Schreiber et al] ilOl Ibl) : [51 lEhnegreen & Elmegreen, ( .2006) : [61 Miralles-Caballero et al. (in prep) 



(z~l-4) have identified massive star- forming galaxies 
with irregular structures where the star formation 
appears to proceed in large clumps, with sizes in 
the 0.1-1.5 kpc range. The most recent are the 
surveys conducted in the GE MS and GOODS fields 
dElmegreen & Elmegreenl20() 6; Elmegreen et al. 20()7p, 
the U l tra Deep Field (UDF: Elmegreen & Elmegreenl 
20051 Elmegreenl boOTt Elmegreen et al.l 120091) and 
the Spectroscopic Ima ging survey in the Near-infrare d 
with SINFONI (SINS: lForster Schreiber et aDl201 Ihh . 
While some of these galaxies show s tructural charac- 
teristics of being involved in mergers (lElmegreen et al. 



20071 : Forster Schreiber et alj 1201 lal). most systems 



appea r as thick turbulent disks (lElmegreen et al. 
20091) . Faint tidal features are difficult to discern at 
higher redshifts because of cosmological dimming. It 
is important therefore to make a direct comparison 
between the star-forming knots in (U)LIRGs and in 
clumpy high-z galaxies in order to understand the 
similarities and differences of the star formation in 
the present Universe and at cosmological distances. 

The clumps in high-z galaxies are less numerous 
and intrinsically brighter and more massive (by a 
factor of 10-1000) than the largest star-formi n g com- 
plexes in local galaxies (lEfremovl 1995 ; 
Elmegreen & Elmegreenv005. Elme green et alJl2007 . 
2009; Forster Schreiber et al. 2011b). The integrated 
luminosity of these clumps usually accounts for about 



10-25% of the total luminosity. The larger masses 
in the high-redshift galaxies are thought to be the 
result of high turbulence in the disks. In fact, velocity 
dispersions of 40 km s"' have been derived from 
spectros copic observations of UDF6462, a z=1.6 
galaxy dBournaud et al.l l2008ai This large-scale star 
formation can also be reproduced by means of numer- 
ical simulations which require velocity dispersions 
several times higher t han in local quiescent galaxies 



(IBournaud et al.l l2007h . Finally, these high velocity 



dispersions have already been observed in clumps 
in local ULIRGs dMonreal-Ibero et alj IIoOtI) . since 
gravitational torques in major mergers can also induce 
star formation in large condensations of gas. Yet, the 
last mention has to be taken with care since strong 
winds and shocks in mergers can also induce large 
gradients in the velocity fields (i.e., high velocity 
dispersions). 



In the standard ACDM cosmology (iSpergel et al. 



20071) the angular size is nearly constant over the red- 



shift range 0.7 < z < 4. At closer redshifts, the angular 
size decreases approximately linearly with increasing 
redshift. A kpc-sized clump corresponds to 3-4 pixels 
in the ACS images at high redshift values. Thus, 
clumps much smaller than this could not be measured 
at high redshifts. On the other hand, if there were 
clumps significantly larger than this in nearby galaxies, 
they would be easily discernible. The lack of spatial 



37 



resolution at high redshifts can make a complex of star 
formation look like a clump since the light from the 
different clusters combines, blurs and some contrast 
with the local background is lost. Whether the high-z 
clumps correspond to a single entity or to complexes of 
star formation, the largest complexes can be identified 
and compared in both local and high-redshift galaxies. 
The high-z galaxy clumps have masses typically be- 
tween 10^ and 10^ M© from z = 1 to 4, decreasing to 
about 10^ Mq at z = 0.1-0.2. The star- forming regions 
in the few interacting clumpy galaxies observed have 
sizes similar to star complexes in local interacting 
systems (i.e., the Antennae or M51), and their mass is 
10-1000 times higher On the other hand, complexes 
in the Hickson Compact Group 31 (HCG 31), stellar 
complexes, which are sensitive to the magnitude of 
disk turbulence, have both sizes and masses more char- 
acteristic of z = 1-2 galaxies (Gallagher et al. 2010l). 

In section l5!2.1l we have seen that knots in ULIRGs 
are intrinsically brighter than knots in less luminous 
systems. This means that for complexes of similar 
sizes and similar population (table |4] shows that the 
colors of the knots in LIRGs and ULIRGs are sim- 
ilar) we are sampling knots in ULIRGs (they could 
also be understood as clumps or complexes of star 
clusters) that can be more massive than in our low 
luminous LIRGs by at least a factor of 4. Although 
the knots in our ULIRGs only account for 0.5-2% of 
the total luminosity (much less than clumps at high 
redshift), we can compare them with the complexes in 
high-redshift galaxies in order to determine if the star 
formation is similar or not to that in clumpy galaxies 
. Akeady at z ^ 0.05-0.1 a significant percentage of 
knots in ULIRGs are more luminous than Mf^i4w=-15 
mag (see Fig. \T5[ . According to the models used in 
this paper, if we assume that these knots are formed 
by a single young (e.g., t~6 Myr) population they 
can be as massive as few lO'* - few 10^ M©, with a 
median value of 1.4x10^ M© . An older population 
(e.g., t~50 Myr) implies a knot with higher mass 
by one order of magnitude. Correction for internal 
extinction makes the population even more massive. 
Thus, these values have to be understood as a lower 
limit. In order to compare the knots in these ULIRGs 
with the clumps at high-redshift galaxies, we applied 
a Gaussian blur to our images in such a way that 
we simulate their appearance at z= 1 (the scale at this 
redshift corresponds to 7.69 kpc/"), using the same 
technique as explained in section |43] 

The knots in the blurred images (hereafter high- 
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Fig. 20. — : Mass vs projected distance for complexes 
in simulated ULIRGs at z=l (in black for a population 
of 6 Myr and in gray for a population of 50 Myr) and 
clum ps in z=2 star forming gala xies (red; values taken 



nps_ 

iJfH 



from lForster Schreiber et alj201 Ibh 

z knots) were measured. Table |7] show the results 
together with the values for clumps in galaxies at 
high redshifts and TDG candidates found in our sam- 
ple of (U)LIRGs. Some of the high-z knots were 
not recognizable in the blurred images because they 
were too faint or too small; a few were blended with 
other knots that were too close or with the nucleus. 
Most of the lowest mass knots and those close to 
the nucleus disappear, although the remaining high- 
z knots range from 2x10^Mq to 7x10^ Mq, with a 
median value of ~ 2.5xlO''M0. The effective radius 
averages about 700 pc and the total size is typically 
> 1 .5 kpc. Taking into account that the derived masses 
correspond to lower limits, they are typically one or 
two orders of magnitude less massive than clumps in 
high-z galaxies, having similar sizes. If the popula- 
tion that dominates the light of the high-z knots were 
older (e.g., 50 Myr), the masses of the high-z knots 
would significantly approach those of clumps in high- 
z galaxies. 

Having set a lower limit for the mass of knots in 
ULIRGs, the mass of kpc-complexes in these systems 
is not likely to lay far from clumps in high-z galaxies 
with a similar size. Therefore, the extreme cases of 
low -redshift star-forming galaxies (ULIRGs) are likely 
to induce large-scale star formation similar to that 
observed in high-z disks (and some few mergers) with 
a high gas fraction. 
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We have also investigated how the projected dis- 
tance of the high-z knots relates to their mass (see 
Fig. |20] While other works do not find any clear 
relation, we measure lower mass high-z knots with 
increasing distance. Although, according to our results 
it is clear that the mass is mainly concentrated at 
small distances in ULIRGs, we can only compare with 
one study of high-z galaxies and we explore larger 
distances by a factor of 10. 

6. Summary and conclusions 

A comprehensive study of bright compact knots 
(most of them apparent associations of clusters) in a 
representative sample of 32 (U)LIRGs as a function 
of infrared luminosity and interaction phase has been 
performed using HST -ACS B- and I-band imaging, 
providing hnear resolutions of 10 to 40 pc, depending 
on the distance of the galaxy. The sample has been 
divided into three infrared luminosity intervals (low, 
LiR < 11.65; intermediate, 11.65 < Lir < 12.0; and 
high, Lir > 12.0) and in 4 interaction phases, based on 
the projected morphology (first approach, pre-merger, 
merger, post-merger). We have extended toward lower 
infrared lumin osities pre vious stu dies focus ed only on 
ULIRGs (e.g. JSurace et al. 199&ISurace et al. 2000) . 
and have detected close to 3000 knots, larger than 
a factor of ten more than in the previous studies of 
ULIRGs. The main conclusions derived from this 
study are the following: 

- The knots span a wide range in magnitudes 



(-20 < Mf 



< -9 and -19.5 < Mm35w < -7.5), 



and color (-1 < Mf4o,5w-Mf^i4w S 5). The median 
values are M/r435vt' =-10.84, M/r8i4iy =-11.96 and 
Mf43sw-Mfsi4w = 1-0, respectively. 



The knots are in general compact, with a median 
Tetf of 32 pc, being 12% unresolved and a few 
very extended, up to 200-400 pc. With these 
sizes, in general (and particularly for galaxies 
located at 100 Mpc or more) the knots constitute 
complexes or aggregates of star clusters. Within the 
resolution limits there is no evidence of a luminosity 
dependence on the size. A slight dependence on the 
interaction phase, in particular in the post-merger 
phase, needs confirmation with larger samples and 
better angular resolution. 



A non-negligible fraction (15%) of knots are blue 
and luminous. Given their color (Mf435vi'-Mf8i4vv <0.5) 
and luminosity (<M/7435vi/>=-l 1.5) they are likely to 
be young (ages of about 5 to 30 Myr), almost free of 
extinction, and with masses similar to and up to one 
order of magnitude higher than the Young Massive 
Clusters observed in other less luminous interacting 
systems. 

Unlike for isolated young star clusters in other sys- 
tems, we find a clear correlation between the knot 
mass and radius, MocR^, similar to that found for 
complexes of star clusters in less luminous inter- 
acting galaxies (e.g., M51 and the Antennae) and 
Galactic and extragalactic giant molecular clouds. 
This relation does not seem to be dependent on the 
infrared luminosity of the system or on the interac- 
tion phase. 

The star formation is characterized by B- and I- 
band luminosity functions (LFs) with slopes close 
to 2, extending therefore the universality of the LP 
measured in interacting galaxies at least for nearby 
systems (i.e., Dl <100 Mpc), regardless of the 
total luminosity (i.e., the strength of the global star 
formation). This result has to be taken with caution 
because the LF of our knots (which can be mainly 
formed of complexes of star clusters) may not reflect 
the same physics as the LF of individual star clusters. 
Nevertheless, there are slight indications that the 
LF evolves as the interaction progresses becoming 
steeper (from about 1 .5 to 2 for the I-band) from first 
approach to merger and post-merger phases. 

Taking into account distance effects, knots in high 
luminosity systems (ULIRGs) are intrinsically more 
luminous (a factor of about 4) than knots in less 
luminous systems. Given the star formation rate in 
ULIRGs is higher than in less luminous systems, 
size-of-sample effects are likely to be the natural 
explanation for this. 

Knots in ULIRGs have both sizes and masses char- 
acteristic of stellar complexes or clumps detected 
in galaxies at high redshifts (z>l), as long as their 
population is about 50 Myr or older. Thus, there is 
evidence that the larger-scale star formation struc- 
tures are reminiscent of those seen during the epoch 
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of morphological galaxy transformations in dense 
environments. 

Knots in systems undergoing the post-merger phase 
are on average redder and more luminous than 
those in systems in earlier phases. Scenarios 
involving either a different spatial distribution of 
the obscuration or a different evolutionary phase 
in which small knots do not survive for long are 
considered. The evolutionary scenario could also 
explain the trend in size toward higher values in 
post-merger systems. A comparison with numerical 
simulations to further investigate these scenarios 
will be presented in a forthcoming paper (Miralles- 
Caballero et al. 201 1, in preparation). 
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